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On the Differentiability 
of Functions of Several Variables 


JAMES SERRIN 


Let w=u(x%)=u(x%,...,%,) be a real-valued locally summable function 
defined on an open region R of the m-dimensional number space #. It is our 


_ purpose here to investigate the differentiability of « when the gradient of the 


distribution associated with w is a measure. By this we mean that there exists 
a Borel measure «, finite-valued on compact subsets of R, such that the formula 


(1 —fugedx=f pda 


' holds for any continuously differentiable function gm with compact support in 


R. Here p, denotes the vector {@q/dx,,..., dg/dx,}, and the measure « is of 
course vector-valued and unrestricted as to sign. For convenience, a function 
satisfying the preceding condition will be called weakly differentiable. 

In order to discuss the differentiability properties of a summable function, 
whose values are in the nature of things somewhat arbitrary on a set of measure 
zero, we make the following definition: A function (x) will be called essentially 
differentiable if there exists an equivalent function v(x), defined on R, which has 
ordinary partial derivatives {0v/@x,, ..., dv/éx,} almost everywhere. 

Suppose now that w is essentially differentiable. If v, and v, are two functions, 
each one equivalent to w and having partial derivatives almost everywhere, 
then obviously the partial derivatives of v, and vy, agree almost everywhere. 
We may therefore define the derivative of « with respect to x,, 1<k<n, to 
be the class of functions equivalent to 0v/0x,, where v is any function equivalent 
to w and having partial derivatives almost everywhere. This derivative will be 
denoted by du/dx,. Although ow/0x, 1s an equivalence class, one naturally 
thinks of it as a function in the same way as elements of L, are thought of as 
functions. Finally, if w is essentially differentiable, we set 


U, = {Ou/Ox,,..., 0UfOxy\. 


The main result of the paper is now the following 


Theorem1. Jf u is weakly differentiable, then u is essentially differentiable. 
Moreover u,—«', where «' is the general derivative of the set function a. 

The proof of the first statement in Theorem 1 occupies Sections 2 and 3 of 
the paper, Section 1 being devoted to a simple but interesting preparatory lemma. 
At the end of Section 3 we digress from the proof of Theorem 1 to obtain an 
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immediate demonstration of a known theorem for the case when u is continuous. 
The remaining part of Theorem 1 is then proved in Sections 4 and 5, following 
which we briefly discuss the relation between weak differentiability and bounded 
variation. The methods of the paper have their partial genesis in the work of 
GoFFMAN, KRICKEBERG, and DE Vito, and our main theorem complements the 
results of these authors!. It may be worth adding that the only background 
necessary for the present work is a knowledge of the elements of the theory of 
integration, and the Lebesgue decomposition theorem. 


Remarks. There are several immediate consequences of Theorem 1 which may 
be conveniently noticed at this point. 

1. By use of a well known fact about distributions [//, pp. 37—38], it follows 
from the first part of Theorem 1 that if T is a distribution whose gradient is a 
measure, then T is defined by a function having partial derivatives almost 
everywhere. 

2. An essentially differentiable function obviously has approximate partial 
derivatives almost everywhere. Therefore, by Theorem1 and a result of 
STEPANOFF [10, p. 300], a weakly differentiable function w has an approximate 
differential almost everywhere. On the other hand, neither w nor any equivalent 
function need have an ordinary differential at any point. The situation with 
respect to the existence of a differential is discussed more fully in the concluding 
sections, 8 and 9, of the paper. 

3. Finally, by virtue of the formula w,=«’, one sees that w, is locally summable 
and transforms as a Cartesian vector. 


1. A preparatory lemma 


The results of this paper are ultimately based on the well known theorem 
of Lebesgue which states that any monotone function of a single real variable 
has a finite derivative almost everywhere. The following lemma generalizes this 
result to functions of several variables. 


Lemma 1. Let f(x) be a measurable function defined on an open interval I 
in 2". Suppose that f(x) is monotone on almost all lines parallel to the x, axis. 
Then the partial derivative of/0x, exists and is finite almost everywhere in I. 


Before proving this result, some comments are in order. Since Lebesgue’s 
theorem for one variable also applies to functions of bounded variation, one 
might conjecture that the hypothesis of Lemma1 could be weakened from 
“monotone on almost all lines’ to simply ‘“‘bounded variation on almost all 
lines”. This change, however, would destroy the lemma: There exist measurable 
functions of uniformly bounded variation on ail lines parallel to the x, axis, 
say, which do mot possess a partial derivative 0f//0x, almost everywhere in J. 
As a matter of interest, such an example is constructed in Section 7. 


? Although our results might seem to follow from Theorems XII and XIV of [3], 
this is not actually the case. In particular, in [3] each component of Ux is treated on 
the basis of a different function v(x), and, more important, the derivatives OUulex, 
which appear in [3] are defined by the special operation lim A u/Ax,, where AWG 9s = SUITS 
cf. footnote 6 of [3]. ee 
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Proof of Lemma. It is clear from the hypothesis that df/0x, exists and 
is finite almost everywhere on almost all lines parallel to the x, axis. It does 
not follow from this, of course, that d//0x, exists almost everywhere in I. We 
must therefore resort to further argument. 


Hor convenience,» the Set of variables (47, °-., 43, %i4.++-, x,) will be 
menotetl sby Vy, and the Set (44 52.., X2150, X%pagn-os, %q) by (fy). Now for 
x=(x,,v)EI we let 


i=" f(a wea lity Sup eer tle = TO, y) (m = integer) 
be an upper right “‘sequential derivate” of 7. It is clear that E* is measurable, 
and, as one easily deduces, finite almost everywhere in J. The same result holds, 
of course, for the lower right derivate E,f, and the left derivates E~} and 
E_f. Now by Fubini’s theorem 


JIB TE. fldx=J Ls |Bt—E, tldmldy =o, 


since the inner integral vanishes for almost all y. It follows at once that both 
limits 
lim bp (%,+1/m, y)—Ht (Fe, y) lim f(*R—- 1/m, y)—f (rp, y) 
m—> oo 1/m : m—> 00 —1/m 


exist and are equal almost everywhere in J. Let E(x) denote the common limit 
where it exists; otherwise put E(x)=co. E(x) is obviously finite-valued almost 
everywhere in J. 

Now let JT denote the set of points «€J such that E(x) is finite and /(x) 
monotone in x,. If x€ T and 1/m+15Ax,<1/m, then 


tat tlmt1,9)—f m9) — At © teattlm, 9) I(x, 9) 
in 1/m a AGG 1/m-+1 : 


/\ 


provided / is monotone increasing at x (when f is decreasing the inequalities 
must be reversed). Since both the right and left sides of the above inequality 
tend to E(x) as moo, it follows that 0//0x, exists and is finite at each point 
of T. But T includes almost all points of J, and the proof is thereby completed. 


2. Proof of Theorem 1. First part 

Consider any fixed open interval J in #" whose closure is contained in R, 
and let k be a fixed integer, 1S. In order to prove the first statement 
of Theorem 1, it will be enough to show that there is a function v(x), equivalent 
to « and independent of k, such that 0v/0x, exists almost everywhere in J. The 
first step in this process is to obtain an estimate of the variation of w along 
lines parallel to the x, axis. 

Let us associate with each positive number / a corresponding function 
K(é) =K(é, h) with the properties: 

i) K is non-negative and smooth for all values of €=(,..., &,). 

ii) K<Const./h" for |E|<h and K=0 for |&|=h. 
ili) f K(&) dé=1. 


24* 
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The € integration in this and subsequent formulas is assumed to be carried out 
over the entire 2-dimensional number space. For x€J (and for h sufficiently 
small), we define the function 


(2) u, = u,(x) =f K(E — x) u(é) dé. 


In view of the properties of K, it is easy to prove that lim , exists and equals 
almost everywhere in 7.2. We set w=lim u, wherever the limit exists, and other- 
wise leave the function w undefined. Clearly w=w almost everywhere. 

It will be convenient at this stage to introduce some simple notation. We 
consider the interval J to be the Cartesian product J J’ of an interval J of 
the variable x, and an interval J’ of the complementary variable y= (41, ..-, %,—1, 
BW Neots Hal Correspondingly, the set of-variables. (x,; +i; %p21, bo %,258 es 
will be denoted by (¢, y) as in the preceding section. Binaliy, for any set B 
in #” we let B, denote the set of points in B whose distance from #” — B is 
greater than h. 

Now by differentiation of (2) and use of (1), there results 


6) jet =— [ Fe Eu de= f KE—2) dey, 


J 


where «, is the k= component of the measure «. An easy application of Fubini’s 
theorem then yields 
J hese oe dx seis) 
For y€ J’ let us define 
V(y) = lim inf if ou 


J h 


Clearly V(y) is measurable on J’, whence by Fatou’s lemma 
Fi Viyydy s fa |a,| 


(integrate first over a compact subset AK of J’, and then let K->/J’). Now let 
p(x) be a real-valued function defined a.e. on J. For y€J’, let Var m be the 
function defined by 


(t, v) | al, he {1/m}. 


[Var ~] (y) = Total variation @ (¢, y) 
ASHI 
with the convention that if m(«) is undefined at a certain set of points, then 
these points are simply to be disregarded in calculating the total variation. 
In view of the definition of V(y), it is clear that the variation of the function 


w= lim ™, satisfies Var w<V(y) for y€ J’. This completes the first step in the 
proof. 


2 Indeed, we have 


Fi Const. "3 
é Sa [ we &)—u(x)| dé, 


[E—x|<h 
and the last expression tends to zero at each point of the Lebesgue set of w. 


| x (*) —u (x) =| JS K(E—~x) [w (8) u(x) 
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3. Continuation 
The function w(x) constructed in the preceding section has the property 


(4) VarwSV(y), fV(y)dy<o. 
Fe 


Moreover, w is equal to u almost everywhere in J, hence is summable. Such 
functions allow a canonical decomposition into the difference of two monotone 
functions, according to the following lemma. 


Lemma 2. Let w(x) be a summable function defined a.e. on an open interval T, 
and there satisfying with respect to some variable x, the condition (4). Then there 
exist two summable functions w* and w~ on I, each monotone increasing in the 
variable x,, such that almost everywhere in I 


(5) w=wt—w-. 


Remark. A representation of the type (5) is not always possible if the equality 
is to hold everywhere in I. Indeed, if it were possible, then by Lemma 1 any 
function satisfying the hypotheses of Lemma 2 would have a partial derivative 
dw/dx, almost everywhere in J. But according to the example cf Section 7, this 
is not the case. 

Proof of Lemma 2. We shall first construct a function w(x), equivalent to 
w(x) but having “‘fewer’’ discontinuities. Let S denote the set of points y¢€ J’ 
such that 

i) w(x,, y) 1s defined for almost all x,€ J, and 

ii) V(y)<o. 

Obviously S includes almost all points of J’. For each x=(x,, y) CJ let us then 
define 

AO 2 ee yeS, 

0 yeaS. 


w(x) is measurable, since it is the limit almost everywhere of the sequence of 
measurable functions {w(x,—1/m, y)}. Now on almost every line y=, the 
functions w and w differ at most at a set of linear measure zero. From this fact, 
and since w and w are both measurable, it follows easily that # =w almost every- 
where in J. The argument above shows not only that w(x) is measurable, but 
also, by a slight variant, that each function @ (x,, y) is measurable in y for fixed x,. 


Now let x,=a be a point such that 
faa, y)dy<o 
y 


(such a point exist, by Fubini’s theorem). Also, for fixed y, let t(x,, y) be the 
indefinite total variation of the function w(x,, y), normalized so that t (a, y)=0. 
We assert that the function 1t(x)=t(x,, y) 1s measurable on I. 

To prove this, denote by P,, the partition of the interval J by m equally 
spaced points t;<f,<-:-<t,,. The approximate variation of w(x;,, y) correspond- 
ing to the partition P, is then 


PACE y) 1 |o (¢;, y) —o(t;1, y)| ‘sign (x; ae a), 
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where the sum is extended over all partition points ¢; which lie in the interval 
between a and x,. Since » is continuous from the left, it is easily seen that 


T (Xp, Y) = lim Tm (X_,¥)- 


Obviously 7, (x,, vy) is a step function in x, for y fixed. More specifically, for %, 
between any two consecutive points of P,,, we have 


T,(X,» ¥) = finite sum of measurable functions of y. 


Thus 1,,(%) =Tp (x, ¥) is a measurable function, and since t(x)= lim T,,(%), so 
also is t (x) measurable. 
We can now define 


w(x) =4ir(x)+o(a)], w(x) = Fir (x) —o(a)]. 


Clearly w* and w~ are measurable, monotone in the variable x,, and w* —w =w 
almost everywhere in J. Furthermore 


|w*|+|w|s {Var w* + $|a(a, y)|} + {Var w +2|a(a, y)|} 
= Varw+|o(a, y)| SV(y) +|o(a, y)|, 


so that wt and w~ are summable. This completes the proof of Lemma 2. 

We now return to the proof of the first statement of Theorem 1. Consider 
the functions (w*), and (w~), defined as in equation (2). Let S’ denote the set 
of points y»€ J’ such that lim(w*),—=w* and lim(w-),—w almost everywhere 
on the line y=y,. Clearly S’ includes almost all points of J’. 

Let us fix attention upon any point x=(x,, y) for which yCS’ We assert 
that if lim (w*), does not exist at x, then lim (w~), does exist at this point. Indeed, 
since (w*), and w* are monotone, if lim (w*), does not exist at x then the relation 


b, Slim inf (w*), < lim sup (w*), S bd, 


holds at x, where 6, and 0b, are the left and right hand limits of wt at x. Thus 
w* is discontinuous at «. But then w~ must be continuous there (recall that 
w* and w~ are the positive and negative variations of a function continuous from 
the left), and hence lim (w-),=w~. 


We are now in a position to define the function v(x) promised at the beginning 
of the proof. In particular, put 


v(x) = lim sup w, (x). 


Clearly v is equivalent to w and independent of k To show that v is differentiable, 
we introduce functions w(x) and w,(x*) by means of 


Pie OI ay Coit lim inf(w-), yes’ 
v(x) y¢S’, 0 vq S'. 


w, and w, are equivalent to w* and w~, respectively, and are monotone on almost 
all lines parallel to the x, axis. Furthermore 


v(x) = w(x) — w(x), x ET, 


Functions of several variables 365 


as follows from the identity u,=(w* —w~),=(w*), — (w~), and the assertion of 
the preceding paragraph. By Lemma 1, w, and wy are differentiable with respect 
to x, at almost all points of J. The same is therefore true of v, and the first part 
of the proof is completed. 

Digression. The preceding proof when applied to a continuous function 
leads to the following important result. 


Theorem 2. If u 1s continuous and weakly differentiable, then u has partial 
derivatives almost everywhere. 

Proof. If wis continuous, then lim w,=w uniformly on any closed subregion 
of R. Thus in the preceding construction v=, and the partial derivative 0u/éx, 
exists almost everywhere. Q.E.D. 


4. Proof of Theorem 1. Second part 


Let v(x) be any function equivalent to w which has partial derivatives almost 
everywhere, and let & be a fixed integer, 1 <<. We must show that 


(6) 0v/0x, =a, almost everywhere in R. 


In the proof, which will be somewhat indirect, we shall maintain the notational 
- conventions adopted in Section 2; furthermore, for simplicity in printing, we 
shall use the tensor notation p, for a partial derivative 0g/0x,. 

| From here on, let F(p) be a non-negative convex function of the single real 
) variable ~, with the property that 


| (7) Fo@+q)s 
_ equivalent to the condition F(p)<|p|-+-const. For any interval I whose closure 
is contained in FR, we assert that 
(8) Je Ee ACRE 
In 

To begin with, for any continuously differentiable function (x) one has the 

relation pat, p+ 
Fidglds) =F (4 f oat nat)s 5 | Feode 
: i 


Xk 


by Jensen’s inequality; we have put 6=A~, for convenience in printing. Inte- 
grating this inequality along the interval J;, where J=J x J’ as in the preceding 


sections, yields 
apt+o 


hoa )dx,< sf fF (p,)dtdx, 


Js *k 


st f fre. ne 


int-s6 


Integrating again over J;, one then obtains 


(9) f F(Ag/Am) dxs J F(p,.) 4 


Ts 
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Setting y=w, in (9), replacing J by J,, and then letting 4->0, we deduce from 
Fatou’s lemma that 


[ F(Au/Ax,) dx < lim inf f F(a, ,) dx. 


I 5 Ih 


Now let w be replaced by v in the left-hand side of this inequality, and allow 6 
to tend to zero. This proves (8). 


5. Completion of the proof 
Here we shall show in contrast to (8) that 


(10) fF(u,,) 44s f F(a,)dx+O(D), 


where @ is a singular measure. 


The measure «(E) can be expressed by means of the Lebesgue decomposition 
theorem in the form 


a(E) =f a'(x)dx+0(Z), 
E 
where #(E) is a singular measure. Accordingly, equation (3) can be written 
Une = JS K(E — x) (6) dé + f K(E — x) dd,. 
Therefore by virtue of (7) and Jensen’s inequality 
F (um, ,) SF(f KE — x) 04 (€) dé) + | [K(E — x) a9, |, 
< f K(E — x) F(oy) d& + f K(E — x) a] 0,|. 
Integrating over J, and applying Fubini’s theorem yields (10), with O=|#%,|. 


We can now prove (6) and thus complete the proof of Theorem 1. By com- 
bining (8) and (10) one obtains 


JF(v,) dx Sf F(a,) dx +O(1). 
i i 


Dividing each side by |Z| and letting J tend regularly to x gives, on account 
of Lebesgue’s differentiation theorem [/0, pp. 118—119], 


(41) F(v ,) F(a) almost everywhere in R. 
Now choose Fs b 
e OF 
F(p) =f 6) = 
De A = Pe Oy 


which certainly satisfies the imposed conditions. Since (pf) is monotone in- 
creasing, it follows from (11) that v ,<«;, almost everywhere. Similarly, setting 
F(p) = (p) = (— 2), we find in the same way that v ,=a; almost everywhere. 
Combining the last two inequalities completes the proof of Theorem 1. 


6. Bounded variation and weak differentiability 
A function u(x) will be called of bounded variation in the sense of Cesari if 
there is an equivalent function w (x) defined a.e. such that, for every open interval — 
I in the interior of the domain of definition of u, and for each variable x4 | 


— — 


Functions of several variables 367 


1<k<n, we have VanwaViy)s » GET. 


for some function V(y) which is summable on the interval J’. The function V(y) 
may of course be different for different directions x,. 


Theorem 3. A function u is weakly differentiable if and only if it is of bounded 
variation in the sense of Cesari?. 


Proof. The “only if” part of this theorem has already been proved in Section 2. 
Suppose then that wis of bounded variation in the sense of Cesari. By Theorem 5.4 
of reference [6], the function w is locally summable. 

Now let C(/) denote the linear space whose members are real-valued continuous 
functions on J, with compact support in J. Similarly, let C’(Z) denote the space 
of continuously differentiable functions with compact support in J. For w€ C’(Z) 
and for any fixed k, 1<k<1, set 


A(@) = UDP, aX. 


Evidently A(qg) is a linear operator on C’(J). Moreover, A is continuous in the 
topology of C(I). To see this, let us observe first that Var w<oo for almost 
all y¢ J’. Thus for almost all y there exists a corresponding measure f, on J 
such that 


Xk 
w(x r=JaBy, and fadip,| V9). 
It follows that 
A(g)=J | Sug nd) ay=S|f pap, ay, 
diet EMS Uaaaf 
and so 
JG) SPs line) Malia, 


proving our assertion. Since A(qg) is continuous in C(/), it may be extended to 
a continuous linear functional on C(Z). But then by the Riesz theorem there 
exists a Borel measure «, such that 


A(g) =— J p day. 


This equality, however, is precisely the condition that « should be weakly dif- 
ferentiable on J. Since J was any interval in R, it follows easily that w is weakly 
differentiable in R. 


Theorem 4. I} u is of bounded variation in the sense of Cesari, then for any 
interval IC R and for any k, 1SkSn, we have 


[V(y)dy=J|aujdx,| dx. 
df I 


Proof. By the preceding theorem, w is weakly differentiable in J, and the 
measure «, even satisfies 


Jalal = V(y) ay. 


3 This result is due to KrRICKEBERG; the present proof, which in part parallels 
that of KrRICKEBERG, is given mainly for the sake of completeness. 
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On the other hand, it is clear that 
Jf d|,| = fae dx=f\|u,|dx, 
I I I 

and the theorem is proved. 


7. A counterexample 

For the purpose of this section we shall suppose that J is the unit square, 
0<x<1, 0<y<1, in the plane of the two real variables x, y. As indicated in 
Section 1, we intend to construct a function g(x, y) on I with the properties 

i) g=O almost everywhere in J, 

ii) Var g=1 on each line y = constant,, 
and nevertheless such that 

iii) @g/@x does not exist at a set of points whose outer measure is one. 


The construction is based on the fact that there exists a set O in I which has 
outer measure one, and yet has only two points in common with each straight 
line parallel to the x-axis (see [/2]). 

With each point Q€O let us associate a sequence of points Q,, which lie on 
the same horizontal line as Q, have rational abscissas ~,,, and tend to Q as mow, 
By omitting certain of the points Q,, and then renumbering the sequence, we 
can suppose furthermore that 


(12) 0< >) |0—Qn| <1 


and that the sets {Q,,} corresponding to the two different points of O on any 
horizontal line are disjoint. The function g(x, y) can now be defined. First, 
at any point (x, y)=Q,, associated with a point Q€O, we set 


g(%, 9) =|0 — On|/4% 


where 4'=2'(Q) denotes the sum of the series (12). Then, at all other points, 
weput ex, y) == 0, 

The resulting function is zero so long as » is irrational, hence property i) is 
assured. Property ii) is likewise immediately clear. Finally, at any point QCO 
we have 


lim inf z 


Ag 


= (0,  lingsap ra = 4/42. 


Thus 0g/0x does not exist at any point of O. Q.E.D. 


Remark. Although Lemma 1 is thus shown to be incorrect under the hypo- 
thesis of bounded variation, a slightly modified version remains true, namely 


Lemma 1’. Let f(x) be a measurable function defined on an open interval I 
in 2". Suppose that f is of bounded variation on almost all lines parallel to the 
x, axis. Then there exists an equivalent function h(x) which has finite partial 
derivatives almost everywhere in I. 


The proof is based on the fact that f(x) can be decomposed into two monotone 
functions as in the proof of Lemma 2. We omit the details. 
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8. The existence of a differential 
We have already remarked in the introduction that a weakly differentiable 
function need not have an ordinary differential at any point. It is the purpose 
of this section to investigate the precise situation with respect to the existence 
of a differential. This can best be done by introducing the concept of strong 
differentiability. A locally summable function w(x) will be said to be strongly 
differentiable if there exists a vector-valued function w(x) such that the formula 


(13) —fugdx=Jpypax 

holds for every continuously differential function w with compact support in R. 
Equivalently, w is strongly differentiable if it is weakly differentiable and the 
associated measure is absolutely continuous. By Theorem 1 a strongly differ- 
entiable function is essentially differentiable and has the derivative u,=y. The 
following result now holds. 

Theorem 5. Suppose that u(x) is strongly differentiable and that u, belongs 
to Ly for some p>n. Then u is continuous and has a differential almost everywhere. 

The conclusion of the theorem refers of course only to a particular member 
of the equivalence class of w. 

For the case of two independent variables, Theorem 5 is essentially due to 
CESARI [2]. The proof we give is a modification and simplification of his. Theo- 
rem 5 is best possible: If w is strongly differentiable and w, belongs to L,,, then 
u need not be continuous, and, even if it is continuous, it need not possess a 
differential at any point. Examples illustrating this behavior are given in the 
following section. 

Proof of Theorem 5. It is well known that if w is strongly differentiable and 
u, isin L,, then there exists a sequence of continuously differentiable functions 
U,(x) such that 
(14) Slum —uldx+>0, ff |ty,~—u,|?dx +0, 

s s 


where S is any compact subregion of R.4 Now let x and y be two points of R, 
and let G denote the sphere whose diameter is the line segment joining x and y. 
If G has compact support in R, then 


(b>), 


where the constant depends only on » and ~. A simple proof of (15) may be 
found in reference [8], pp. 523—524. 

From (14) and (15) one easily infers that the sequence {w,,(x)} is uniformly 
bounded and equicontinuous in any compact subregion of R. Consequently there 
exists a subsequence which converges to a continuous function v (x) in R. Obvious- 
ly v=u almost everywhere, and the first part of Theorem 5 is thereby proved. 

We may assume from here on that w is a continuous function, and that 
U»—>~u uniformly on any compact subregion of Rk. Letting m-—>oo in inequality 

5) then yields 


fi 1 
(16) |«(x) — w(y)| < Const. |x — y|*-"/? (J [ee]? dx) is 
( 


(1 5) | Mn (x) — Un (y)| < Const. | x _— gy ame (J | ohn xl” dx) 


Hence u(x) is Hélder continuous in R, a result due to Morrev.) 


4 See reference [4]. We can in fact take u,,=u, with h=1/m. 
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Now let E denote the Lebesgue set of w, relative to the class L,.? We shall 
show that « has a differential at each point x, in E, and this will prove the 
theorem. Indeed, consider the function 


w(x) =u(%) —A-x%, 


where A denotes the constant vector u,(x,). Evidently w(x) is strongly differ- 
entiable and w,=u,—A. Therefore, as one sees by comparing with (16), 


1/p 
Jw (x) — w(x9)| Const. [x — gf" f J. — Al?) 


it being assumed that x is near enough to x) so that G is in R. Since % is in E 
the integral on the right is o(|x— x9|"), whence 


| w (x) — w(%»)| = 0 (|x — xl). 
This is precisely the condition that w have a differential at x). Q.E.D. 


Remark. Since a Lipschitz-continuous function is strongly differentiable and 
its derivative belongs to L,,, we obtain as a corollary of Theorem 5 that every 
Lipschitz continuous function has a differential at almost all points of its domain 
of definition (Rademacher’s theorem). 


9. Two examples 
Consider a set of points Q,, k=1, 2, 3,..., dense on some interval I. Let a, 
be a sequence of positive numbers such that >») a,<oo, and define f,(x)= 
log log (1-++-1/r,), where 7, denotes the distance from x to Q,. Finally set 
(17) u(x) = 2 ag fy (%)- 


Since /, 1s summable, it is clear from Lebesgue’s theorem that the partial sums 
u,,(x) of the series (17) converge in L, to u(x); this shows, by the way, that u 
is finite almost everywhere. Moreover, since 


lsl=(Slhsl"@x) <0, 


we see that the sequence {w,,,} is bounded, and hence weakly compact, in L,. 
We let y denote a weak limit of this sequence. Then from 
i) Uy Px ax =— J P Un x ax, 


it follows by letting m-—>+oo through an appropriate subsequence that 


Jug,dx=—fowdx. 
Thus w is strongly differentiable on J, and its derivative is in L,. But wu is 
infinite on a dense set, and so is everywhere discontinuous. 


. We shall next construct a continuous function w which has strong derivatives 
in L,,, but which nevertheless is nowhere differentiable. (In the paper quoted 
earlier, CESARI proved the existence of such a function in the case nN=2, but 


: See [14, p. 65]. The work here refers only to the one-dimensional case, but 
carries Over to m dimensions in the obvious way. 
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_ he did not give an explicit example.) The method is an elaboration and extension 
of one due to STEPANOFF. 


| The function to be constructed will again have the basic form w=» ap te 
_ but will be more complicated in nature. We set 


Uj R/2*, Nee SUMEATI a ae (t) = log £5 “ [log ee 


where e=e, is fixed by the condition log e=—2*. Also let 0, be a positive 
number less than 27°", to be fixed later. In order to define f, (x) we partition 
the interval J, which may be assumed to be a cube of unit length on each of its 
sides, into N;’ subcubes each having side length 2~”. Let Xz) be a generic symbol 
for the center of any one of these subcubes. We then set 


I, (x) =A 4 FE %4)|/Ox) 


if |x — x%)| So, for some center x), and otherwise put f,(*)=0. (In geometrical 
language, the surface /, consists of a family of ‘“‘cones’” of height one and base 
0;,, centered at the lattice points x); in two dimensions this is an Indian mystic’s 
bed of nails.) 


Since 0Xf,<1, it is clear that the series >’ a, f, converges uniformly to a 
continuous function u. Moreover w is strongly differentiable. Indeed, the partial 
sums w,, satisfy 


mM 


[b= [3 af 


by Minkowski’s inequality. But 


<> ale 


: \1/n 
It | an N, lo | = Const. N, ( i | co’ (t) |” jr dt) 
0 
S Const. N, | log eels, 


as one easily sees. Consequently ||2,,|| << Const.) 2*°-”!", so that {w,,,} forms 
a uniformly bounded sequence in L,,. By our earlier argument, then, w is strongly 
differentiable with w, belonging to L,,. 

It remains to show that « cannot have a differential at any point. Consider 
first a point € where all the functions /, vanish. This point lies no farther than 
\/n|2*** from some vertex %@). Hence 


u (a(p)) —U (&) > Ck fr (*(x)) = 2k 
[@aesl... Vue y 


Since the right-hand side can be arbitrarily large for appropriate choice of xq), 
it is evident that w cannot have a differential at €. Next, if only a finite number 
of functions 7, are non-zero at &, an easy modification of the above argument 
again shows that w cannot have a differential at &. 

Finally, suppose that f, (€) +=0 for an infinite number of values of k.. Consider 
a fixed k for which /,(€)=-0, and form the difference quotient with respect to 
the nearest of the vertices xj), or else (if |%)—&|<%o,) with respect to an 
appropriate point yg) where |x.) —¥q@)|=@x- Then, provided the radii @, have 
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been chosen sufficiently small, one sees that. the dominant contribution to the 
difference quotient is due to the function a, f/,. But the slope of this function 
becomes arbitrarily large as k tends to infinity; hence again w cannot have a 
differential at €£. Thus ~ does not have a differential at any point in J. 


Note. After this paper was completed, the author learned that C. GOFFMAN 
and R. E. Hucues had independently established that a function of bounded 
variation in the sense of Cesari is essentially differentiable. Their proof (which 
will be published elsewhere) and the work of Section 3 of this paper are closely 


similar. 


The above research was supported by the United States Air Force Office of 
Scientific Research under Contract No. AF 49(638)-262. 
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Il Teorema del Massimo Modulo per I’ Equazione 


dell’ Elastostatica T ridimensionale 


Dedicata ad ANTONIO SIGNORINI nel Suo 70° compleanno 


GAETANO FICHERA 


Sia A un campo (insieme aperto) dello spazio cartesiano delle tre variabili 
reali x1, %,, 3; considerata in esso l’equazione vettoriale omogenea dell’equilibrio 
elastico 


(1) A,v +k grad divv=0 


2 
(4: = >» na , v vettore a tre componenti reali; k costante reale maggiore di — 1] ; 
eat 


Ci proponiamo, sotto opportune ipotesi per la frontiera X' di A, di dimostrare 
Vesistenza e pervenire al calcolo di una costante H dipendente unicamente da 
A e da k, tale che per ogni soluzione della (1) di classe due in A e continua in 
A=AvS riesca: 

(2) max |v (x)| S H max |v(x)|. 


/ 


Risultati come quello testé enunciato sono classicamente noti per le equazioni 
del secondo ordine di tipo ellittico in una funzione scalare e si conoscono gli 
analoghi teoremi per l’equazione parabolica del tipo del calore. Recentemente 
si sono ottenute delle estensioni alle piu generali equazioni di tipo misto ellittico- 
parabolico [2], [3]. Mrranpa ha considerato alcuni anni addietro |’equazione 
biarmonica ed ha conseguito per essa un teorema di massimo modulo [10}, e, 
pitt recentemente, ha esteso il suo risultato ad un’equazione ellittica di ordine 
2m in due variabili [77]. AGMon [7] ha quindi considerato il caso di equazioni 
ellittiche in pit variabili. 

Occorre pero dire che in queste ricerche non é stato preso in considerazione 
uno degli aspetti di maggior rilievo del problema, specie dal punto di vista delle 
applicazioni, cioé il calcolo effettivo delle costanti, analoghe alla H che inter- 
viene nella (1). Tale questione, che non si pone per le equazioni ellittiche o 
paraboliche del secondo ordine in una funzione scalare, riuscendo allora H=1, 
si presenta assai complessa negli altri casi. Per quanto riguarda la (1) é stato 
solo osservato da G. Pérya [13] che per k=-0 non puo essere H=1. Inoltre é 
stato conseguito [4] il calcolo di una costante H¢ tale che per ogni soluzione 
continua della (1) sussista la seguente maggiorazione 


v(x), 


max |v (x)| S Hp max 


essendo C un qualsiasi compatto contenuto in A. Tuttavia la H¢ ottenuta é 
tale che per CA riesce im He =+~. 
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1. Ipotesi, premesse e tappe della dimostrazione 


Supporremo che A sia un campo limitato e connesso la cui frontiera 2’ sia 
costituita da »-+1 superficie regolari chiuse: X (contorno esterno) e 2;,..., 2, 
(contorni interni) a due a due disgiunte, ciascuna di classe Ci. Con cid s’intende 
che ogni X, é dotata di piano tangente in ogni punto e che esistono tre numeri 
positivi ¢,, ¢, # (0<AS1) tali che, fissato comunque un punto x di ', accade 
che: 7) l’intersezione Q,, (x) di X' con la sfera di centro x e raggio 4 € un insieme 
connesso; 77) assunti come asse x; la normale interna a 2’ in x e come piano % %2 
il piano tangente, l’insieme ,,(x%) @ una porzione di superficie di equazione 
X3—=@(%, X») con P(X, %2) definita in un dominio T,? del piano x, x2 (proiezione 
ortogonale di Q,,(x) sul piano tangente); 77) la funzione p(x, %,) ha derivate 
prime verificanti una condizione di Hélder uniforme in 7, con coefficiente c ed 
esponente h. 


Se B é un insieme dello spazio cartesiano a tre dimensioni, dotato di punti 
interni, con C’”(B) indicheremo la classe delle funzioni reali di classe m in B, 
cioé dotate di tutte le derivate parziali fino all’ordine m, ciascuna delle quali 
coincide con una funzione continua in B. Dicendo che una funzione vettoriale 
appartiene a C”(B) intendiamo che ciascuna delle sue componenti é una funzione 
ci Gia) 

Per ogni vettore u(x)={u,(x), u2(x), us(x)} della classe C1 (A), detto n il 
versore normale interno a 2 e dies con d/én loperazione di derivazione 
secondo tale versore normale, porremo: 


n+ (14 7) oe4 ae 
Y Bal bl On. * Boar 


L(u) =|(k ~ 5*;)diva 


k 
ceo: (n A rot u) 
[(a@ 6) denota il prodotto esterno dei due vettori a e b|; k é wna costante reale 
maggtore dt —1. 
Porremo anche: 
E(u) = Anu +k grad div uw. 


Diciamo WU la classe dei vettori reali a tre componenti appartenenti a 
Cl(A)OC?(A), nulli su Ye tali che E(u) appartenga a C°(A). Denoteremo con 
(A) [-¥Y(2)] la classe dei vettori w(x) tali che |w(x)| sia pseudo-limitato in 
A (su 2), cioé tale che, tolti i punti di un insieme di misura nulla in A (su %, 
rispetto alla misura superficiale), la |w(x)| sia limitata. Indicheremo con p.e.s 

A 
| w (x ) |[p-¢ e.s |w(x)|] lo pseudo-estremo superiore di w(x), cioée l’estremo inferiore 
dei numeri Q tali che |w(x)|<@Q quasi ovunque in A (su 2). 


Consideriamo per ogni w€ ll l’equazione 
(1.1) f pL(u) do= f + E(u) dx 
sy A 


(do = misura elemento d’area su , OX=04%, a2, d Xs). 
Supponiamo che sia assegnata in ¥°(X) e si ricerca in L(A) una w 
incognita verificante le (1.4) per ogni well. Tale y esiste comunque si assegni 


1 
Re amine intendiamo la chiusura di un insieme aperto. 
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g se e solo se esiste H>0 tale che per ogni w€Ul riesca 


| (4.2) J\Loo|dosH f |EWlaxs 


_ Se sussiste la (1.2), esiste almeno una y tale che 


(1.3) — -—pesly| SHpes|gl. 


Se proviamo che la soluzione yw della (1.1) é unica, assumendo quale @ la traccia 
su 2’ di una soluzione v della (1) appartenente a C!(A)-C2(A), deve necessaria- 
mente essere (per l’identita di Green) y=—v e quindi la (1.3) altro non é che 
la (2) dell’introduzione. Per provare tale diseguaglianza per le soluzioni apparte- 
nenti a C® (A)nC?(A), l'idea che si presenta spontaneamente é quella di appros- 
simare uniformemente in A una soluzione siffatta con una successione di soluzioni 
appartenenti a Cl (A)OC2(A). Occorrerebbe a tale proposito un teorema di 
densita rispetto alla convergenza uniforme delle soluzioni di classe uno in A 
nello spazio delle soluzioni della (1) ivi continue. Questo teorema sussiste se 
si riesce a provare che ogni soluzione della (1) appartenente a C°(A)C2(A), 
nulla su 2 é identicamente nulla. Tale circostanza sara da noi dimostrata. 
D’altra parte, per conseguirla, si provera che sussiste per ogni soluzione v della 
(1) continua in A ed ogni wEU, Videntita 


(1.4) fuvL(u) do=— f vE(u) dx 

>) A 
talché la (2) si potra dedurre direttamente dalla (1.2) tramite la (1.4), come nel 
caso delle soluzioni di classe uno, senza far ricorso al teorema di densita anzidetto. 


E opportuno indicare i vari punti dei quali constera la dimostrazione della 
(2). Avvertiamo che con uw indicheremo costantemente un elemento di U. 


I) Esistenza e calcolo di due costanti positive H, e H, tali che: 
(1.5) max | L(u)| S Hy max | E(w)| + He ( flu|taa)). 
II) Dimostrazione della diseguaglianza : 
(1.6) ( J\ultax) <H, (. f|E@Pax) 


con calcolo esplicito per la costante positiva H, e, quindi, esistenza e calcolo di 
H,>0 tale che 


(1.7) max | L(u)| = Hy max| E(u). 


III) Esistenza e calcolo di H,>0 tale che: 


max |u| < H,max | E(w)|. 
A A 


7 


1.8) 
La (1.8) costituisce un teorema di massimo modulo per le soluzioni dell’equazione 
E(u) =f non omogenea con la condizione omogenea al contorno: w=0 su 2%. 


2 Cfr. [2]. 
Arch. Rational Mech. Anal., Vol. 7 25 
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IV) Dimostrazione della completezza nello spazio C°(A) (munito della con- 
vergenza uniforme) del sistema di tutti i vettori {E(w)} e conseguente deduzione 


(per dualita dalla (1.8)) della diseguaghanza 
(1.9) HN Ra Bes Gia 

V) Esistenza e calcolo di due costanti positive Hg e Hy tali che 
(1.10) (L@ldes Heda pias |ulae, 


e conseguente deduzione della (1.2) da (1.10) e (1.9). 
VI) Dimostrazione della (1.4) per ogni soluzione della (1) appartenente a 


C°(A)C?(A) e conseguente deduzione della (2). 


2. Dimostrazione del punto I 
Poniamo a,=—4(2n+3)(2n+5)(2u+7) essendo mn un intero positivo 


arbitrario, e 
t 


yO=a,2e—1yt, FQ) = f <r — ape ar. 


2 
1 
Riesce: 
£(0) =0; EO) ==" F’"(0) =0, FOU) =0 (t= OA evan 56 
aii 
t FY t dt os (2n+3) (2n a i 5) (2m+7) = b 
(2.1) / eu 8 (m+ 1) (m+-2) (m+ 3) 1 
Lin) = 47, 
Sige 
(2.2) 1 ae [FO OJ Oprah yc. 4, Fhe 


Sia ¢ un fissato numero positivo e, in corrispondenza ad esso, poniamo: ins 
aa) per 0SrSe, f(r)=0 perr>e. Si verifica che / (0) =0, f’(0) =4, /’’"(0) =0, 
f?(e)=0 @=0, 4, 2,3, 4; 5) e quindi f€C5(0, +00). Riesce inoltre: 


(2.3) Ml<tR, Sipr®|rar=s,. 


Siano %= (x, %2, X3)€ ¥=(y,, Ve, Vg) due punti dello spazio. Posto r= |x—-y|, 
siha perr>0: 


ef (7) = f''(r) (%;— Vi) (4a—n) + f’(r) 


Ox; OXp, 
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Poniamo: e2 } 
Be ON po apt ek 
} y af Ox; Ox; ji Mk 2(1+h) ‘ 


sel yY oe i rt 
Ej = 0) Aa — te paras ean 


Indichiamo con S(x, y) la matrice 3 x3 il cui elemento generico é la funzione 
5; (%, y) = E% (vy). Essa € una matrice soluzione fondamentale per l’equazione (1) 
ed € nota col nome di matrice di Somigliana®. Si ha pi (x, y) = si (x, y) + Bi (x, 5) 
con $;(%, y) =O(1). 

Sia p'(x, y) il vettore di componenti #/ (x, vy), £3(x, v), p(x, y) e sia L, 
(4=1, 2, 3) la h-esima componente dell’operatore differenziale al contorno L sopra 
introdotto. Si ha: 


ee ee (ee es 9 dat, 
et ae, 2 pee iat RACE ey Fag Gena he Be n 
= ay ae Bg ‘| k é i | 
fp an* k+2 on* y| ! 
aa (4%; — i) (4,—Vn) tt 0 1 
(2.4) raat —f so pe i d / (x; Vi) (x), Vn) on® 72 | 
1 OY (%i;—Yi) (%n—Vn) y 0.1 
a eS i 
1 7? (4%3—9;) ny, + (*n,—Vn) Nj # 
a) Vy (3 
Si ha d’altra parte per ogni punto & di A: 
— 4 u,(€) = f L(u) 2’ (E, y) do, +f E(u) p (5 y)dy — fuk [p(s »y)]ay. 
Posto E'[w): Ei [w,] [w=(w,, Ws, Ws)], si ha E’[p*(x, y)]=4 0; AGH (7), sicché, 


posto g' = (6; MARS. O AR) iA ry basitelec’ 
— 42u;(§) = JEW PE, y) doy + SE) b (Ey)dy — af usilé y)d 


Sia x un punto di 2’, n* la normale (interna) a 2 in x, € un punto di A su n’. 
Applicando ad ambo i membri della precedente relazione l’operatore L relativo 
alla normale n* e facendo tendere € a x su n* si trova (cfr. [5] teor. XIII) per 
ogni xE€ 2: 

(2.5) —22L(u) =f L(u) *{Li(G)} do, + [EW * {Lie )}dy—3 [ux Li@)jay 


Ne viene: . 
2x |L(w)| Sm2x|L(w)| J Li) ||4o, + 


(2.6) is max | E(w)| f Li (pay +4 allie) [Lie ay. 


Grice, aul 

S Scrivendo L* si intende che l’operatore L opera su p'(¥, y) come funzione di 7; 
con n* si intende lanormale interna a Y in 7; 71, %, 23 sono le componenti del versore 

normale interno a 2. . - 
> Con il simbolo w* {aj} indichiamo il vettore di componenti aju;, aju;, aj uj. 
2 2 
Porremo inoltre || @}|| = D(a)?’ - 
i,j ; 


By: 
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Supponiamo, cosa ben lecita, che il numero 94, in precedenza introdotto quale 
costante geometrica legata ad A, sia minore di 1. Assumiamo e<gy4. Sia 2, (x) 
l’'intersezione di 2X con la sfera chiusa J,(x) di centro x e raggio «. Poiché 1 
vettori #’ hanno, come funzioni di y, supporto contenuto in J,(x), si ha: 


{Mew ldo, = Fae )|ae,, JIE \ay = f \ZE@ bay. 


Per la (2.4) si ottiene: 


k yap 
fitie Pde, V3 I aa én* Ast R42 On® 7 #Ssap 
Ss 
|R| tr Pee vt 4) tr laa ie: 2) 1 | i 
ss k+2 [ll | o — 7 fe / én* | on® doy 
Q2¢(*) 
[al aaa V3, || ii a Ser 
- ee) ey Ido, =(F aE a Ue ine 
2¢(%) 2¢(x) 
2V3+3|R| ee at ah it ne 
pal tle e Hi sn ir" 
2¢(x) Hae 
Si ha: 
| mas as | (%1—1) Put (%2 _ Px. — (%3— Vs) (1-4 Pz,+ Px) 3 


<= |(4- V1) Px, + (%2— Vo) Pag — 9 (%1» Xe) + (M1, Vo)| 


=17|(%1— V1) Pe, + (%2 — Va) Pre — Pra [41> Xe +8 (Yo — ¥9)] (%2 — Yo) — 
— Px,[%1 + O' (9, — %), Yo] (4 — 44) | 
Sr la — y,| cr +72|%— y2| er? S2c7*. 


Sia @, # un sistema di coordinate polari nel piano tangente a Y in x, avente 
poloan %, Si ha: 


doy=(1+ G5, +95,)' dy do S V1 +207" od dd < (1+ |2¢ 0") pdoad 


e pertanto, tenendo presente la prima delle (2.3): 


J[IEE@ do, < 4 e(1 +20) (3 + [AI )3 fot do+ 
2 0 


2 k+2 


= 213+318l | Ras . 
smn femoan): 
0 0 


= |k| £, . h+. k : mt 
+2)2me 7 a] e detan det fi (e)|e de 
0 


0 


[eset aTone 


+ 4n0(1+ 2c) 215s rah (= ,R+4 | “F) +20 A fire )| odo. 
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Ne segue, ricordando la seconda delle (2.2) 


7 U k . wer 
(2.7) pie (p))|do, < M, ®, & + 2218 fli (e)lede 
0 
con 
wre = \ st IR Pte 
™ M,=4me{[(1 +V20) (7 ! ol ba tol ati ! 
ay = ONG eset hes 1 
t (ict V2e) C er male 


Posto A,(%)=I,(x%)OA per la (2.4) riesce: 


fswrters (S85) furl 


2343 |A| ” 
+ 21S lacs 
A,(x) 


da questa, tenendo presente che 


ay + 


ee 


iv [itp 


oe 


o 


oe fir 
A 


e(%) 


Fea A; [ -4y Sane, [ pray SAne, 
si trae: Ae(x) A¢(x) 
V3, 21A| ; 
pies pay < ~i | sd | 
Quindi: 
(2.9) fe [Zi (e')|4y S Mp, 
con 
(2.10) iM, = eee) pou stS Dy € 
3 
Si ha: 
(of fe /2|k ' 
nies (13 +1 +24) eo 
1 
- 2|k Ps V2|kl\ & 
= (3+ lal + BEY Morte ferenad ss (B+ + as) 
0 
donde 
£ fiwjizeeilays4(_ f 1et@ Pay) ( fi ray) 
A Ag(x) 
Te 3 
< V3z (V3 +|2| + ae lias i fires) 
Ne viene 
(21) Ma( fluleay) 2 2S lol ILHe Ney 
con 


(2.12) My = V3 (V3 +14| +2) Rent. 


380 GAETANO FICHERA: 


Da (2.6), (2.7), (2.9), (2.44) si trae 


2a max |L(u)| <M, ®; a” max | L (uw) 


aa 


& 
ee 


ode)mgs L(w)| + Mz max |E(w)| +34 flwtay). 
—- 4 wy } 


Tenendo presente la seconda delle (2.3) e la (2.1), si pud assumere m tale che 


lel eet bee Let ie: 
(2.14) 1— sob > S(t Ep My ay] ani EB, 


Una volta fissato 1, nel modo indicato, si scelga ¢ tale che 


(2.45) eS \5 
cid che implica e<o,4. Si ha allora: 


2n(1 = maar b,| yi lae Se 


(2.13) si trae: 


Poiché pud assumersi ®, = 


nN |? 


jaz (1 — aan b,) —M, |a,| el max | L(x) 


2+hk 
E(u) | [ 


A 


Si deduce allora la (4.5) con 


M, M. 
(2:46) Hy iF] 3 : Re iF 3 a 
il ear bn) —My [an |e! 2a(1— en by | —My | ay, |e” 


dove M,, M,, M, sono date da (2.8), (2.10), (2.42). 


3. Dimostrazione del punto II 

Proveremo la (1.6) come caso particolare di una formola di maggiorazione 
relativa a sistemi di equazioni pit generali di quello dell’elastostatica ®. 

In questo paragrafo A denotera un campo regolare dello spazio cartesiano Sy 
(g=1), cioé tale che 0A =@4A, essendo la frontiera 0A costituita da porzioni i 
ipersuperficie regolari in guisa tale che sussistano le formole di Gauss-Green 
che mutano integrali estesi ad A in integrali ipersuperficiali estesi a GA. Sia m 


un intero positivo e supponiamo che le funzioni reali ai (x) (BAe eons 
t,J=1, 4», m) del punto *=(%,,..., %,)di A appartengano a C2(A), le bi; (x) 
a C1(A) e le c;;(x) a C°(A). Poniamo: 
Dae a) hp paee o? 1 eee AR TD) bP D 
he Ox ) hike a OXp Ox, , 1 a a; j hk aE Mie? +: Chis 
EX, = D, aj} D, — D, bi, + ¢;;. 


° C/r. [4|, dove la maggiorazione in questione trovasi data per la prima volta. 
Cfy. anche [12]. 
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Facciamo le seguenti ipotesi: 


hk kh hk git Ye ae AES 
) aj; =4i;, af = Aji» Or = U;53 


1 

2) la forma quadratica a} } (x) A} i} & semidefinita positiva per ogni x € A; 

3) la forma quadratica c;,(x) 4,4, & semidefinita negativa per ogni «€ A; 

4) esiste una funzione v €C?(A) non positiva in A e tale che la forma quadra- 
lica E¥(v) A; A; & definita positiva per ogni x€ A. 


Per ogni coppia di funzioni v e w di C2(A) si ha: 
Gut) J [oF ,,(w) —wEK@)]dx=— f [o(at} Dy) m —w(at} Dyv-m —b,;2)] do, 


essendo n, la k-esima componente del versore normale interno a 0A e 5;, 
(2), — —D, att) M,. Si assuma nella (3.1) come v la funzione indicata nella eee 
4) e w=u;u,;, essendo u=(u,...,4,,) un vettore appartenente a C1 (A) C2(A) 
nullo su 0A e tale che E;,(u,;) appartenga a C° (A). Poiché riesce: 


: E(u; uj) =u; E;;(u,) + u; E; (uj) + 2ai* Dj, U; Dy Uj — C55; U;, 
ers yhalcs E;,(u;) +4, £;;(u,)]) dx + 2foai; Dj, u,;- Dy, uj, 4% — 
ie SP i! NEE Nas — 0. 
Per l’ipotesi 1) la (3.2) puo scriversi: 


(3.3) 2fvu, E,j(u) dx+2 J var; D,,4;-D, uj 4 %— ae jit a Sry ER (e) da 
A A ; 


€ per le1);'2);-4),, 4) riesce: 


fue u;-D,ujdx — foe, u,ujdx So. 
: A 


Dalla (3.3) si trae allora: 
fu;,u,ER(v) dx 2 J vu; E;;(u;)dx. 
A y 


Indicato con Q il determinante della matrice {E*(v)} e con Q,; il complemento 


™m 
algebrico dell’elemento E7;(v) di tale matrice, riesce: E¥j(v) A; A; 2 Qm> Y O77 J. 
osteo max (12... a<, faq.) segue: Joie 


ie Toi ee ge | IulPaascons [09 Oi ax 
; 
=) it vu, E;;(u;) dx < 2max |v| [|m |E(u)|dx, 
; A 
avendo indicato con u il vettore (2, ..., %m) € con E(u) il vettore (E,,(u,), ..., 
E,,;(u,)). Abbiamo quindi: 


(3.4) ( f jul2dx) < 2mmax (Q* Q>) max f JEw[eae), 
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Sia ora g=m e E(u)= A,u-+-k grad div u, (k>—1). Detto d il raggio di una 
sfera che contiene A nel suo interno e di centro x®, si assuma leer dors %o |? — a2. 
E facile verificare che sono soddisfatte le 1), 2), 3), 4). Poiché riesce: E;;(v) = 
E* (v) =26;(m+hk), Q* Q4=[2(m+h)]+, dalla (3.4) si trae: 


( flees) eral flee (w|tax) 


E quindi provata la (1.6) con 
} ge 
(3.5) Cr guy; ? 
essendo d il raggio di una qualsiasi sfera che contiene il corpo elastico. 
Dalla (1.5) e dalla (1.6) si deduce la (1.7) con 
(3.6) H,=H,+ H, H,(vol A)i. 
4. Dimostrazione del punto III 


Per ogni x€A si ha 


— 4% u(x) al L{u(y)] * {s;(x, y)} doy +J E[u(y)] * {sj(*, y)}dy. 
Si trae: 
4n|u(2)| Smax |L(a| { [si(x, 9) do, + max (9 f Isi(x, »)]dy. 


Le funzioni di x: 
1 i 1 i 
aq [li Nlde,, Ze [Ist vlay 
s A 

sono continue in S,; dettine H; ed H, i massimi in A, sussiste la diseguaglianza : 

max | (x (x) = A, max |L(u) u)| + He max | E(w u)|. 
Da essa, tramite la (1.7), si deduce la (1.8) con 
(4.1) =H, A, ay. 


Per maggiorare H, e H, si osservi che 


Isis = foi ay) £— mad io) 
1,3 ea 
— 1 , z 3 k 2 k 1 
S15 eal > + [> (ee 9 ape] = BEML 
Detto % il punto di massimo su & della funzione if ny / Si-hax 
doy % ‘fares S2n0 Wi pet eee ee 28 
re vy = a pai [9 os ‘ 
204%) 2 — Qo (2) 
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Riesce d’altra parte | --—- <2 02, essendo 6 il diametro di A, per modo che 


Ze a 


possiamo Baie a H, ed H, al modo seguente: 


R k x 
(4.2) Hy ao Ses el | (3 2 V14 2c — 2), 
V3(k+2)+|k| 0 
(43) Hes Gadhia 


5. Dimostrazione del punto IV 


Sia w(x) = (w,(%), we(x), ws (x)) un vettore appartenente a §°(A)?. In par- 
ticolare, potremmo assumere come w(x) un vettore le cui componenti sono 
polinomi. Consideriamo nella classe dei vettori continui su »' l’equazione inte- 
grale nell’incognita : 


(54) 27-0 (% + J Ply) ) « {L? (s*(x, y))} doy TN ) « {L¥[s*(x, y)]}dy =0. 


Poiché riesce L,(s'(x, y))=O(|~—y|"~*) per x, yEX, la (5.1) & una equazione 
di Fredholm (di II specie). Se la soluzione di essa esiste, appartiene a §°(2’)8. 
L’equazione omogenea trasposta di (5.1) é: 


(5.2) 2a (x) + fry) *{L} (s(x, y))} do, =0. 
= 
Si consideri, in corrispondenza della soluzione t di (5.2), la funzione: 


x) = J o(y)* {L? (s*(x, y))} do,; 


essa appartiene a $'(A)%. Poiché z é in A soluzione della equazione dell’elasto- 
statica e dato che si annulla su 2 essa é identicamente nulla in A 2°, 
Poiché riesce: 
Jx(x) do, f(y) * {Lj (s (x, y) ))}dy= sf 10 (4) 2 dx=0, 


resta provato che la (5.1) ammette soluzione. 
Se m é una soluzione di (5.1), poniamo: 


(53) u(x) =—*- [ oly) #{S)(x, y)}do,—  [ wly) {5a )} ay. 


La funzione vettoriale u(x) gode delle seguenti proprieta: 1) é continua in 
tutto lo spazio; 2) detto CA il complementare di A, u appartiene a ! (A) eda 
§1(C.A); 3) E(u)=w in A; 4) E(u)=0 in CA e, come funzione di §1(C A), u 
verifica la condizione al contorno L(u)=0 su 2; 5) w=0 su 4; 6) wha le derivate 
prime con modulo di quadrato sommabile in C A; 7) é infinitesima all’infinito. 


7 Dicendo che w(%) appartiene a ”(B) (m=O), intendiamo che esso é un vettore 
di C”(B) ele sue derivate m-esime verificano una condizione di Holder uniforme in B. 

8 Cfr. [7], p.170. La dimostrazione ivi contenuta é relativa al caso k=0, tuttavia 
essa sussiste praticamente inalterata per ogni k> —1. 

9 Cfy. [7], p. 212. Sussiste anche in questo caso l’osservazione della nota’. 

10 Cfy. [5], teor. II. 
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Da 5), 6) e 7) segue che w é nulla all’esterno di 2," e, se »>0, € costante 


all’interno di ogni 2X; (i>0). Se é y=0, resta provato che well. Se é€ v>0, 
esistono 3 autosoluzioni linearmente indipendenti della (5.2); esse sono costituite 


dai vettori 7”"(x) (m=1,...,v; 7=1, 2,3) aventi la s-esima componente ont 
cosi definita: Oh Pa Pepe 
Tj (% 
; =0 7eL'— 2, 
Siano g(x), ..., g®”(x) autosoluzioni linearmente indipendenti dell’equazione 
omogenea seve a (5.1) e sia 4 °(x) una ben fissata soluzione di questa equa- 


zione. Poniamo (x)= ¢(x) + 3 g(x), essendo le c,, costanti. Sia w(x) 


definita da (5.3) e u°(x) dall’ sais di (5.3) ottenuta sostituendo m con @g?®. 
La funzione “(x)=u(x)—wu®(x) é nulla all’esterno di 2) e costante all’interno 
di ogni X’, ({>0). Non possono scegliersi le costanti c,, in guisa tale che #(x) 
sia nulla all’interno di tutte le 2’. Infatti, nel caso contrario, sarebbe 

30 


=J( Zeng" (y)) # {5(a, 9)} de, =0 
3p 
e quindi >) c,, p(y) =0 su 2. Ne viene, scelto x” interno a 2,, che la matrice: 
m=1 
{ fo" Oi" »)} do, 


ha determinante non nullo. Pertanto esiste una (sola) 3v-pla di costanti c,, 
tali ae 


s Cm i; p(y Nes {5 (x7, y)} do, 
m=1 2 
= — f(y) *{si(x", ybdy — f g9(y) «{si(x", y)} do, igen 4g Doras 
A Pay 
In corrispondenza a tale scelta delle costanti é w(x”")=0 (n=1,...,¥) e quindi 


u(x)=0 in C(A). Pertanto é wEUl anche per y>0. Dalla proprieta 3) di uw e 
attesa l’arbitrarieta di w segue la proprieta di completezza asserita nel punto IV. 
E opportuno notare che se » é di classe C?, essendo allora la u(x) teste 
costruita appartenente alla classe C2(A), la proprieta di completezza in discorso 
€ verificata dai vettori nulli su 2 ed appartenenti a C2(A). Tale classe, che 
ovviamente é contenuta in U, sara da noi in seguito indicata con Uy. 


Dopo quanto si é provato, pud affermarsi il sussistere della (1.9), ottenuta 
per dualita dalla (1.8)1%. 


6. Dimostrazione del punto V 
Dalla (2.5) si trae: 
Aon oe (y))| doy S La (8 (x, y))||do, + 
(6.1) i 
(JIE HO) dy {Li (P(x, 9) | do, + ES le(o)| ay SIL (e'( »)) doy. 


One ian teor. Ot: 
HO Olle 
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Con procedimento del tutto analogo a quello seguito nel §2, da (6.1) si trae: 


ee |dos|N, |a,| e+ 2218p, | [leo] do+ 


(6.2) 
+N, f Boal ay-+H, f leonlay. 
A A 
ove 
A 4) 4 3 (2 (24-02) A] 20343 /k[ 74 7 154) 
edad dnc }5 2+h pa) ae 2 -EEe ie 2 ronal 


[2 ie 
A) New} 1003 +03414) 18] pTAL 8, 
(6.4) N, 2(2-£h) ag Lhe age ee do,+ 7 he (area 2’) 


; N= V2 || .—3 x [ doy 
(6.5) N,=3(V3 + al + Sal) Re max foe 
Si fissi l’intero m tale che, oltre le (2.14), sia anche verificata la diseguaglianza: 
2) 
IN GE cle 2a(1 Raae by). 


I] numero positivo ¢ sia tale che, oltre a soddisfare la (2.15), soddisfi anche 


la seguente condizione: 
safe AL gy 
—— Qa: 


Dalla (6.2) si trae allora la (4.10) con 


eh ee ae) Cee nee eS 
2a(i——l », |= Niger oF damaiss | 5, — Nap] et 
24k 2+k ” 
Si noti che se, scelto comunque un punto 4° in A e una semiretta di origine 
in x°, questa sega »' in un numero finito di continui disgiunti che non supera 
mai il numero fisso yu, riesce allora 


(6.6) He=- 


Cleo 
= 47 
sup |: . eas | 22 AGE IL 


wed én’ 
Dalla (1.10) e dalla (1.9) segue la (1.2) con 
(6.7) H Hee His: 


Resta cosi provata la (2) per le soluzioni della (1) appartenenti a Ct (A) C2(A). 
1213.6), (4-delG.7)suitae: 


H =H, +H, H;H,+H,H,H,H,(vol A)!+ H,H, 


provvedendo al calcolo esplicito delle costanti a secondo membro le (2.16), (3.5), 


13 Chr. [9]. 
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La costante H si pud far quindi dipendere, oltre che dal parametro k, dalle 
seguenti costanti geometriche del campo A: ¢, h, 04, vol A, d, 0, bm. 


7. Dimostrazione del punto VI 

Supponiamo che 2 verifichi l’ipotesi, pit particolare, di essere di classe Cj. 
Sia A(é) un versore definito in ogni punto & di 2, le cui componenti siano funzioni 
di classe 2 su Xe che formi col versore normale a X interno ad A un angolo < 32. 
L’esistenza di un versore quale /(£) é ovvia. Esiste un @)>0 tale che per 0< a 0 
il luogo XY, descritto dal punto x=&+ @A(é) quando &€ descrive 2, é la frontiera 
di un campo A,, la cui chiusura é contenuta in A. Si ha, detto w un arbitrario 
vettore di C2(A) e v una soluzione della (1) appartenente a C° (A)oC2(A): 


(7.1) —fvE(u)dx+J {uL,(v) —vL,(w}do,=0, 

Ap =o 
avendo indicato con L, l’operatore L calcolato su X,. Per la continuita di v 
in A riesce: 


(7,2) lim fvE(u)dx=fvE(u)dx; lim fuL,(w)do,=fvL(u)do 
00 4, A g>0 3 


Diciamo # un arbitrario vettore di classe 2 su 2’e sia u tale da appartenere 
a C2(A) e da esser cosi definito in A—Aj,,: u[&+oA(&)|=u(é). Indichiamo 
con Y, lo spazio dei vettori (a tre componenti) di classe due su 2’ munito della 
ordinaria topologia della convergenza uniforme celle funzioni e delle derivate 
prime e seconde di queste (rispetto ai parametri locali su 2’). Sia ZF lo spazio 
duale di Z,. Poniamo per 0<o0< 40, 


I, = f wl, (0) do,= f #L, 0) “e do; 


Dy yy 
“= 


“@ 
riesce T,€ DF. Qualunque sia %#, esiste, per le (7.1), (7.2), lim<T,, #>. Ne segue, 
o—0 . 

per il ben noto principio di Banach-Steinhaus della uniforme limitatezza}, 
Vesistenza di un elemento T€ QF tale che bay hee ,4)=<T,m. Segue allora 
dalla (7.1): 

JvE(u) dx + fuL(u)do=<T, uy 

A i 


e, quindi, per w€ Up, il sussistere della (1.4). Dalla proprieta di completezza 
stabilita per Uy nel §5, si deduce, per dualita dalla (1.2), scritta per w€ Up, il 
sussistere della (2) per le soluzioni di (1) appertenenti a C° (4)>C?(A). 
L’ipotesi, ammessa al principio di questo §7, che & sia di classe C?, si potrebbe 
rimuovere osservando che H dipende dalle costanti di Hélder degli elementi 
differenziali del primo ordine di 2, epperd, supposto XY soltanto di classe C?, 
basterebbe approssimare uniformemente XY con una successione »,, di frontiere 
di domini A,, interni ad A, tali che i coefficienti di Hélder c, (relativi ad un 
fissato esponente h) degli alent del primo ordine su ,, si mantengono equi- 


limitati, riuscendo ogni %, di classe C?. Ma di cid ci heaneitauee lasciare la 
cura al lettore desideroso di tale etencoss, 


14 Cfr. (8). 
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(Pervenuta alla vedazione lo 11 gennaio 1961) 


On the Solutions of a Class 
of Equations Occurring in Continuum Mechanics, 


with Application to the Stokes Paradox 


I-DEE CHANG & ROBERT FINN 


1. Introduction 


In this paper we apply methods developed by FINN in [/} to study asymptotic 
properties of solutions of equations 


(1.1 a) Ave b- Vib a 0. D— conse. 
(1.4b) KA = Cp, ¢ = const}, 
for a vector field 26 (x) ={w,,..., w,} and scalar field #(a) which are defined 


in a neighborhood & of infinity in m-dimensional Euclidian space. The notation 
in (1.1) is the usual one of vector analysis. By choosing particular values of 


the constant c and the constant vector b, we obtain from (1.1) the following 
equations of fluid and solid mechanics (see, e.g., [2, 3): 


C—O}, b=0 Stokes equations of hydrodynamics 


Mi G==@ Oseen equations of hydrodynamics 
191) (D==30 Equations of linear elasticity. 


We assume in most of this paper, as is the case in the problems of mechanics 
from which the equations are derived, that c2=0. However, in some contex‘s, 
notably the basic Theorems 1 and 4, it is sufficient to assume c>—1. 

The main result of this study consists in showing that under a surprisingly 
weak growth condition at infinity, the solution admits a representation in terms 
of a single surface integral over an inner bounding surface ». From this re- 
presentation we conclude that the solution tends to a limit (t#), fo) at infinity, 
and the rate of decay of (w#, #) to this limit is controlled by the asymptotic 
properties of the fundamental solution tensor associated with (1.1). We de- 
termine this tensor explicitly in §3. In §4 we apply our result to a discussion 
of solutions for which @(a) vanishes on XY. These solutions are of importance 
physically as steady fluid flows or as deformations exterior to a rigid body. 
We obtain in particular a new clarification of the Stokes Paradox of Hydro- 
dynamics, which states essentially that in two dimensions a solution of the 
(time-independent) Stokes equations in & which vanishes on is necessarily 
unbounded. Our result is sharp in all dimensions and shows that the singular 
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behavior of the solution is intimately connected with the forces exerted on in 
the motion. We obtain as special cases all earlier results on this and on related 
problems which are known to us [cf. 3, p. 361, 4, 5, 6,7]; the demonstrations 
are In our opinion simpler and more natural than those previously employed. 

The methods of this paper are not restricted to equations of the form (1.1), 
but are applicable to any elliptic system for which the fundamental solution 
tensor exhibits the qualitative properties needed in our demonstrations!. Equa- 
tions with variable coefficients can also be treated, and it seems likely that a 
general description of asymptotic behavior can be constructed along lines sug- 
gested by N. G. Meyers [8]. The particular choice of the equations (1.1) has 
been based on the consideration that they are general enough to indicate the 
scope of the methods, but not so general as to cause the ideas to be obscured 
by technical detail. 


We remark that in the range of c considered the system (1.1) is elliptic (the 


characteristic determinant is (1+c) (}é?>0 for c>—1, Y&> 0} but not 
i=1 

strongly elliptic, hence not of a type for which boundary value problems have 

been studied in recent literature. The properly set boundary value problem 

requires prescribing w(a) on the boundary, but not the pressure (a), this 

function being determined by the solution®. Correspondingly, in the material 


1 It would suffice, for example, for the fundamental solution tensor to have the 
property that for any prescribed positive integer k, all derivatives of sufficiently 
high order (depending on k) of each component of this tensor behave at infinity 
as O(y—*). Professor L. HGRMANDER has communicated to us a proof of the follow- 
ing equivalent condition, valid for any system with constant coefficients: Let P(D), 


where D=(—7D,,..., —1D,), 1=imaginary unit, denote the determinant of the 
differential operators occurring in the system, and let, respectively, §=(&, ..., &,), 
¢=(¢,,..-,¢,) be vectors in real or complex Euclidean n-space. Then the fundamental 


solution tensor has the above stated property if and only if both the conditions 
2) erie) fore a0; 


b) ImG>co if $-co onthesurface P(C)=0, 
are satisfied by the polynomial P. The representation theorems of this paper, and 
the corresponding asymptotic estimates of the solutions, can be extended without 
essential change to all systems which satisfy the conditions of HORMANDER. 


Such systems include, for example, the linearized magnetohydrodynamic equations 


Di Ot nA eat 
L | eo| Wo 
V-w=0 


for which the essential features of Theorem 1 are easily seen to apply. In this case, 
there is no Stokes Paradox, since the fundamental solution tensor vanishes at in- 
finity for all n22. 

2 The existence of a solution corresponding to such boundary data has to our 

> 
knowledge been proved in the literature only in the special case c=b=0 (Stokes 
equations). It seems likely that, starting with the solution of Opgvisr [9] for this 
case, the existence in the general case can be demonstrated by the continuity method, 
using estimates of the type studied in [0] and in [JI]. A direct demonstration 
ag . 

by the method of surface potentials in the case c=0, b arbitrary, has been given by 
J. SANDERS in a work to appear shortly. 
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of this paper no assumption is made on the behavior of the pressure, and the 
regularity of the pressure at infinity is proved as a consequence of growth 


. om 
assumptions on 2 (a) alone. 


2. Notation and Definitions; the Integral Representation 

The following notation will be used. We choose the origin O of our coordinate 
system to be a point interior to 2. The set of all points exterior to 2'is denoted 
by &. Let Xp denote a sphere of radius R centered at O and enclosing 2’ We 
denote by % the annular region bounded between 2 and 2. The region 
exterior to X', is denoted by &p. 

Let 6 (a) =[w;(x)] and T(a)=[T;,;(a)] be, respectively, vector and tensor 
fields defined in &. We define 


We write 7;;=O(R*) if there exists a constant A so that |T;,(a#)|< Ala|* as 
a—>co; we write 7;,=0(R*) if |w|~*|7;;(#)|—+0 as aco. Similar definitions 
apply to vector and higher-order tensor fields. 

Corresponding to a given vector field &(a#) and scalar f(a), we define the 
second order tensor field 


eS a Ou; Ou; 

(24) Ti = {Ty = — (1-9 95, + (SH 4 2) 
where pe 
i= |) Grek 
Once a: 


We also define the linear differential operators Lu and M@ by the expressions 
Li=AU—Vp—b-Vi, 
Mt=A¢’—Vq+d-Ve. 


For any vector fields &# and & and corresponding scalars #, q satisfying the 
conditions 


(2.2) V-u=—cp, V-b8=—cq 


in a region V with sufficiently smooth boundary S, we have always the identity 


(23) f@-Lé—&.M8)aV=f [B-TU—W.TE— G-¥) (0- BI as 
V s 
where by & - T& is understood the expression 


v- Tu= v,{T U},,n,, 


n = {n,} is the unit exterior directed normal on S, and summation is extended 
over repeated indices. 
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Let x (@ — y) = {yi}, Y (a — y) = {y;} denote an nxn matrix and an n-vector, 
respectively, which become singular at w=y in such a way that 


din=lim f [Ty +x (b-i)] ds 


where 2, denotes the surface of a sphere of radius 7 about # as center and Tx 
is formed by interpreting the components of W as pressures. We obtain from 
(2.3) the integral representation 


(2.4) t(x) =f (y-Lu—U-My)aV — f [y- TU—&-Ty—(-y) (O-*] AS. 
V Sy) 


In (2.4), # appears as a parameter, the differentiation and integration being 
performed with respect to y. 
Let now w(x) and x(a —y) be vector and tensor fields which satisfy (2.2) 
and also the equations 
Tw) =05 “UM (y)'=0 


for each column vector in the matrix (y;;), with y as independent variable. 
We have then the integral representation for the solutions of (1.1), 


(2.5) w (x) =— f [y-Twe—w-Tx-— (w-x) (b-n)] Sy 
s 
A similar reasoning yields the representation for the scalar f(x), 
(2.6) p(a) =— f [b- Te —wW-TH— (W-) (b-”)] dS. 
s 


In forming TW we have introduced a “‘pressure’’, p*=b : v(-). 


3. The Fundamental Solution 
The fundamental solution (x,y) is determined from the equations 


(3.1) oe ee 

(3.1) At = +B - V4; = 54: 5(Y — 2). 
One puts (cf. [72, p. 31]) 

(3.2a) = — in (A+b-V)®, 

(3.2b) tri = Oust A +ob-V]O— 3 


Equation (3.1a) is then identically satisfied, while equation (3.1b) reduces to 


(3.3) (A+b-P) [1+ 4+cb-V]G=d(y—2). 

If b+ 0; let b=|b|, (= ue —- . A solution of (3.3) can be determined in 
the form 

(3.4) Dly—a) => i, [1+ c) ®,— D,] dr 
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where ®, and @, are solutions of the equations 


3.5) (A+b-V)®,=dy—a); [1+ 4+cb-V] O,=d(y —2). 
Lira); 


n— 1> 

1 Ba bY Me hers) 

Oe ,=— salaer) Keoal a] 4 ‘ 
Lost ea ea) 
Gopi Ss Te ERB Lemp Sleipewy : ; 


a 
while if c=0, the expression for ®, is replaced by 


(3.6c) @, == logr, n=2, 


@, = 7G) ; NS, 


2? (n—2) "2 


(3.6) 


where y=|y—a| and K,_, is a modified Bessel function of the second kind 
Bh me 

(cf. [12, p. 36]); if b=0, equation (3.3) reduces to 

(3.7) (1+¢) 44G=d(y—2), 


and we may take 


(3.7b) =a = : N= 2 ane 4, 
8x3 (2-2)! (140) 
and 
r(}—2| wats 
(3.7¢) De oe. ; nm =odd and »=even> 4. 
16%? (1-+¢) 


When these expressions are substituted into (3.4) and (3.2), the fundamental 
solution tensor (7, %) is obtained explicitly. From (3.6) and (3.7) we may compute 
formally all asymptotic properties of this tensor. However, for the principal 
purposes of this paper we need only the following crude estimates as 7->oo, 
easily obtained from the form of the expression?: 


a) b =0: 
(3.8) Ix] =O(logr), |¥| =O(r>), n=2, 
G.8b) =O), =00-, nda, 


(3.8c) |VYy)=OG7-*->), |VMw) =o"), kE1, nZ2 


where V denotes a k'" iterate of the gradient operator. 


* Somewhat more detailed estimates in the case c—0O are given in [11] and [73]. 
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b) b+ 0: 
(3.9a) XO PEO) a= 2, 
1+k F 2+k 
(3.9b) |V | =O1r I, \POe|=olr =), ID 
(3.9¢) =o ® ), el =olr 2), 23, 


4. Statement of Results 


In this section we state the main theorems of our paper. We prove these 
theorems in §6. Because of the qualitative difference in asymptotic behavior 


which occurs, we discuss separately the case b=0. It is the singular behavior 
of the solutions in this case which gives rise to the Stokes Paradox. The case 


b+0, and some theorems valid for every choice of b, are introduced later. 


Case 1. b= 0. 


Theorem 1 (Representation theorem). Let w(a) be a vector field defined in 
an n-dimensional neighborhood & of infinity and which satisfies (1.1) for some 
scalar field p(x). Suppose that as r->oo, either 

i) w (a) =o(7) 
or 
ii) f |wl?7—*?) dV = 0 (log R) 
IR 


OY 
il) i, | V w|?7-" dV =e o (log R) . 
SLR 


Then there exist vector and scalar constants Wy) and py such that throughout & 


(4.1 a) W (av) = Wy = (x: Twé—w-Ty)ds, 
(4.1 b) p(x) = po a [P- Tw —w-Tylds. 
The fields W (a) and p (x) adnut the asymptotic expressions 
(4.2a) W (x) = Wy + x (x) -E+O(r-"), 
(4.2b) P(e) = by +(x) +00) 
where 

C= — firwds 

55 


Citiee Le 


Theorem 1 is best possible, in the sense that “‘o 
in the hypotheses. We remark that the vector © can be interpreted physically 
as the force exerted on the contour 2’ in the motion or displacement. 

26* 


cannot be replaced by “O” 
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A particular consequence of Theorem 1 is that 


= 
with an improved behavior at infinity in case Z=0. This observation bears 
on our next theorem, in which we introduce a boundary condition on 2. 


Theorem 2. Under the hypotheses of Theorem 1, if w(x)=0 on &, then 3-0 
unless W(x) =0in €. If n=}, then 


(4.3) T- ib) = 2 f (defew)2dV +c(1 —o) f p2aV 
y é 


6 


Ae Ow; Hy. os 
D\ OMS Aaa 

The physical content of (4.3) is that in any steady fluid motion past 2’ with 
adherence on the boundary, the work done at the boundary per unit time equals 
the rate at which energy is converted into heat in the motion. Thus, in three 
or more dimensions no energy is introduced into the flow at infinity. This is 


where def vw denotes the matrix, - 


not the case in two dimensions when b=0. In the interpretation of v6 (a) as 
an elastic deformation, (4.3) asserts that the work done at the boundary due 
to a constant deformation — 2, on XY equals the potential energy of deformation 
in the medium. A particular consequence of (4.3) is that in any fluid flow past 
a rigid body 2, or in any elastic deformation of a medium which surrounds 2’ 
and is stationary at infinity, the body necessarily experiences a “‘drag’’ force 
in the direction of 26). 


Theorem 3. Let w(x) be a solution of (1.1) in &. If for some wo, 


then Fy =() 


From Theorems 2 and 3 we obtain in particular the Stokes Paradox of Hydro- 
dynamics. For if (a) represents a two-dimensional fluid motion exterior to 3, 


and if w#(a) is bounded, then necessarily =0. But if the fluid adheres to Pa 


then 2 (a) =0 on , hence z =- 0, a contradiction. If c=0 our result is a strength- 
ened form of the results in [3, 4, 5, 6]; if c> 0, we obtain a result which includes 
the extension to this case by Durrin & Nott [7]. 


Case 2. b+ (0) 


Theorem 4, Under the hypotheses of Theorem, there exist vector and scalar 
constants W, and py such that throughout €, 


(4.4) (a) == Wy = J [x: Tw —w-Ty— (w-y) (b-7i) dS, 


(4.4b) p(t) = by — Sips Te —w6- Teh — (6-) (b-)] dS. 
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The fields w(x) and p(a) admit the asymptotic expressions 


(4.5) 6 (x) = 1, + Y-M+O(r), el 

(4.5 b) Ww (a) = wy 4 xm or 7, n> 2, 
: = aptetel 

(4.5) Pe) =Po+p me +o(r =), 

where 


M—= — f [Tb —w(b-v)] ds 


and 


> 


me — — (i [Tw = b (wu nN) — Ww (b : n) | as. 


Theorem 5. Under the hypotheses of Theorem 1, if wW(a)=0 on &, then T+0 
unless Ww (x) =0in &. Also, for all n=2, 


(4.6) ¥- i) =2f (def w)2dV +c(1 —c) f praV. 
é om 


The following theorems are valid for every choice of b. The first of them 
includes the work of Finn & Nott [6] and of Durrin & Nott [7] on the same 
question. 

Theorem 6 (Uniqueness theorem). Let w(a) be a solution of (1.1) in & for 
some scalar field p(x). Let w(x)=0 on X and suppose that as > either 

i) w(x) +0 
or 

ii) f|w|?7-" dV =o(log R). 

IR 


Then W (a) =0 in @. 
With any fluid motion for which (a) +2, at infinity there is associated 
the kinetic energy of disturbance from the uniform flow, E= f |w#— w,|2dV. 
é 


In any potential fluid motion, this quantity is always finite and can be used 
as the basis for a variational approach to the study of such motions. The situation 
is, however, quite different in the case of solutions of (1.1), as can be seen from 
the following theorem?: 


Theorem 7. Let n=2 or 3, and let w(a) be a solution of (1.1) which satisfies 
the hypotheses of Theorem’ and for which w(a)=0 on X, w(a)==0 in & Then 
J |e (x) — Wo|2dV =o for every choice of Wy. 

é 


5. Preliminary Lemmas 


Lemma 1. Let u(x) be any Cartesian component of a vector freld Ww (a) which 
satisfies (1.1) in & for some scalar field p (a). Suppose that for all sufficiently large R, 


(5.4) f wr-dV < K(R) 
I R41 
4 A theorem of this type has been proved also for the (non-linear) Navier-Stokes 
hydrodynamical equations, in two dimensions by W. WoLIBNER [/4], and for the 
three-dimensional case by Finn [15]. Important extensions of WoLIBNER’s results 
to time-dependent motions of compressible fluids in two and in three dimensions 
have been given by A. Krzywickx1 [16]. 
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for some real s. Then for any non-negative integer k, 
|V*) « (a)| =O [r? K4(r)] 


where V*) denotes the k* iterative of the gradient operator. 


Proof (cf. [17]). One finds from (1.1) the following equation governing u (a): 


L,u={i+dgA—cb-V}{A—b-V}u=0. 
If one defines the adjoint operator L, by 
L,={(it+oA+cb-V}{A+b-7}, 
one obtains the following identity for any scalar functions and y in a region #: 
(5.2) p Ly —ypLl,g=divW 


where 


=(1+¢) [pV(4y) —ApVo + Ap Vy — yV(4q)] — 
— (1+ 2c) [pV(6-Vy) +yV(b- Vo) —bVp-Vy)] + 
teb[pb-Vy—yb-V@). 


Let P be an arbitrary point of large magnitude in &. Construct a unit sphere 
S with P as center and let Q be an arbitrary point of S. Let €(Q) be an infinitely 
differentiable function defined in S such that 


O=15 ter lo alah 
Let us now put p=u and g=C®@ in (5.2); w and @ being the solutions of 
Lu=0,. £,0=6(P—). 
Integrating (5.2) over the unit sphere S, we find 
(5.3) 4D) = =a (2) Le[S(Q) P(P— Q)] dV 


where # is the annular region }<|P—Q|<3. We can differentiate the right- 
hand side of (5.3) with respect to P, thereby obtaining a representation for the 
derivatives of wu ee 


VO u(P) = i Q) L2[o(Q) Vp? ®(P — Q)] dV. 
Thus, letting 7 denote Artis to the origin of a reference system, 
[V™ u(P)|?s f wr-saV fr ri DAltho D(P— Q)|?dVo 
af af 
by Schwarz’ Inequality. Thus for a suitable constant C, 
|V") w(P)|2?S C(rp +.1)5 f wr-sadV 


SF R+1 
=O[r K(7)], q.e.d. 
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Lemma 2. Let P(x) be a polynomial of degree m in &, and let s be a real number 


such that 
f P?(@) 1-8 dV =o (log R). 
IR 


Then P(x) =0 if sSn, and m<34(s—n) if s>n. 


Proof. We have 
P(@) = a,,(Q) 7" +--+ +a, 


where the coefficients depend only on the position 2 on a unit sphere. We obtain 
JPA ale dV SH or OP dr de Ol ye 1 Yar dO 
ISR IR IR 
=e Soh at i EL dy re O(1) ae at aaa 

for some positive constant k?. 

Suppose s<n. If then m=0, we shall have f P?(@)7r-‘dV> k® log R+0 (1), 

IR 
a contradiction. Thus in this case P(w)=0. If on the other hand s>~z, then 
R 

in order that f7?”-s+"-1d7y=0(logR) we must have 2m<s—n, q.e.d. 


Lemma 3. Let w(at) be a continuously differentiable vector-valued function on 
and exterior to a unit n-sphere X, and let w(a)=0 on X. Then 


f |el?r-@* dV < f|Vewbl2r* av. 
Ir FIR 


Proof. On an arbitrary ray through the center of 2, we have 
R > 
Hf ab |2 7-42) "1 dy = — |w/? 


1 


R 


R 
+ | (w-w,) r2dr 
! 
us She 
<|/ f\elrsar Vii |e,|2774 dr 
1 1 


R R 
ee trne ate t ars hl Vee |tr treater, 
1 1 


so that 


and the result then follows by an integration over concentric spheres. 


6. Proofs of Theorems 


Proof of Theorem 1. By (2.5) a (a) has the representation 26 (ax) = 26, (a) + 6,(x) 
where 


(6.1) wy (w) = — [Ix TH —-Ty]4S, 
(6.2) w, (x) =— f[x-Tw—w-Ty]ds. 
=R 
We show first that #, (x) =y-F+0(7-”), We have, in fact, 
W, (x) = —X (a) - [TU (y) ASy—J {x (#9) —x @)] TH Y)— Y) Tx (H—Y)} aS, 
Pay z 


=x (x) -E+0(7-”) 
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since by the estimates (3.8¢) 
[TY] =O (A= yw tow yon 


and 
x (a — y) — x(@)| =l¥- Vx) 
where 
E=2— ay, Og a S15 
so that 


lx(@ — y) —x(x)| =O”) 


again by the estimates (3.8). 
To estimate 26, (a), we observe first that by Lemmas 1 and 3, 


17 26] =O (r°?) 


under any of the hypotheses of Theorem1 for any s>n-+ 2. (In applying 
Lemma 3, we may clearly assume that w(«#) vanishes in a suitable sphere.) 
From the equation (1.1) we then find |Vp|=|4e|=O(r°) for any such s, so 
that | p(a)|=O(rS*?)). Now 


(6.3) w,(@) = — f [x-Tw—w-Tx| dS 


for any sufficiently large R. Since the left-hand side in (6.3) is independent 
of R, so is the right-hand side. Similarly, we find 


1" g(x) —— f (Py. Tab — a. TPM IAs 
=R 


for any integer k, the integral being independent of R. But by the above esti- 
mates and the estimates (3.8c), we see that for any ees , the integral over 


Xp tends to zero as R->oo. It follows that for any such k, V“) w&, (a) =0, from 
which we find that #, (a) is a polynomial of degree at most k+ 1, w(x) =P,, (x), 
mM=k+1. Thus, 6 (x) = vb, (a) + 0, (wv) = P,, (x) + 0(r) by the estimates (3.8a, b). 
If w(a)=o(r), we conclude immediately P,,(#) =const.=%#,, from which the 
stated representations of (a) follow. If 


i 


IR 
we observe (since |P,, (a)|2> cv? in a semi-infinite cone with vertex at the origin) 
that necessarily 


— 


w 


2y—("+2) dV — 0 (log R), 


JP, Biikre er) eV Sto (leew) « 
FR 


V 2— . => . . 
By Lemma 2, we have m< "te —" =1, ice. P,, (x) =const.=w,. Finally, if 


Y 


z |Vw|?7 "dV =o (log), the result follows from the above discussion by use 
rR 


of Lemma 3. The stated representations for f(a) are obtained similarly. 


i Eeot of Theorem 2. Let us suppose =0. Then by Theorem 1, 2 (a) = 
Wy +O0(r-"), b(t) =po+O(r"), and (by the same proof) Vwé(a)=O(r—"), as 
y—>co. Formal integration by parts, applied to the solution 26 (a) — a,, establishes 
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the relation 


T-Wy)=0=2 f (def W)2dV + c(1 —c) f p2dV + $ ( — &,)- TWAS. 
Fr Fr ZR 

Because of the above estimates, the outer surface integral vanishes in the limit, 

and we find 


(6.4) 2 f (def ®)2dV +c(1 —c) f p24V =0. 
é é 


If c<3, we verify easily, using (6.4), that defw#&=0 in &. Hence the motion 
is a pure rotation. But #=0 on 2, hence 06 (#) =0 in &. This result, however, is 
true for every non-negative c. To see this, we observe that the divergence 
theorem and the above asymptotic properties of w(a) imply (summation on 
repeated indices) 


since w,=0 on 2. Adding this relation to (6.4), we find, using (1.1b), 
[ Veb)2dV +c f p2dV =0, 
& 


and we conclude again #=0 in @. 

If n=}, the relation (4.3) can be obtained by the method used to prove (6.4), 
the outer surface integral vanishing in the limit because of the asymptotic 
properties of the solution. For »=2 the volume integrals which occur are nec- 
essarily infinite when +0, as can be shown by using Theorem 1 and the ex- 
plicit knowledge of x (x, y). 

Proof of Theorem 3. By Theorem 1, it is sufficient to observe that if T+ O, 
then 


= O (log 7) , ee 
OW); N > 23 
This property is easily checked, using the explicit knowledge of x (a, y). 
Proof of Theorems 4 and 5. These theorems are proved by the methods 
used for Theorems 1 and 2, the differing estimates arising because of the different 
properties of the fundamental solution tensor in the case b==0. We omit details. 


Proof of Theorem 6. As in the proof of Theorem 2, we use the identity 


f wav + cf Pav =e — p(w-n)| dS. 
IR In ZR 
According to our assumptions, the estimates of Theorem 1 or of Theorem 4 must 
hold with #,=0. It follows that the outer surface integral vanishes in the 
limit. Hence Vae=0 in @, from which #=const. But w=0 on &, therefore 
w (x) =0 in @, ¢@.e.d. 
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Proof of Theorem 7. We first observe that by Theorem 1 or Theorem 4 
(w=0 on ,) 
46 (a0) = y + x (a) -E + W, (w) 


where 26, (a) has finite square integral over & except in the case n=2, b=0. 
(The proof of this fact requires estimates on x (x, y) which are slightly more 
delicate than those indicated in §3, but is not difficult.) Putting the exceptional 
case aside for the moment, we see that Neen w,|?dV is finite or infinite 


according as i Ix: |? dV is finite or AAitte pe the components of = 


by {f,}, the fueeeea becomes a q We adopt 
spherical coordinates and observe that because of the ee ficenerd se of 
y;;, the integral over a sphere centered at the origin of any term for which 7 +k 
must vanish. Thus, we need only consider an integral of a sum of squares. Each 


of these terms has, however, an infinite volume integral unless ==0, as is easily 


verified from the explicit expression for x. If S=0 on the other hand, we con- 
clude from Theorem 2 or Theorem 5 that # (a) =0 in &. 


Te remains to discuss the exceptional case n= 2, b=O. In this case we find 


that if £-0, the ratio of the integrals, { IX: .Z\2dV and ef 2dV, tends to 
Fn 


uw, 


FR 
infinity with R. Thus we conclude again that Ey |e |?dV =co unless == 0, but 
again T=0 implies by Theorem 2 that v6 (a 0 in 6, q.e.d. 


This investigation was supported by the Office of Naval Research. 
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On the Ouasi-Static Theory of Viscoelasticity 


EDWARD L. REISS 


Communicated by E. STERNBERG 


1. Introduction 


In general, viscoelastic deformations are time dependent. However, if the 
inertial forces resulting from the deformation are “‘small’’, it is customary [/—3] 
to neglect them in the equations of motion. The resulting differential equations 
are those of the “quasi-static” theory which may supply an approximation to 
the dynamic theory. For many viscoelastic materials, neglecting the inertia 
terms removes the highest order derivative with respect to time from the differ- 
ential equations. Therefore, one or more of the initial conditions of the dynamic 
theory may not be satisfied by solutions of the quasi-static differential equations. 
Hence there is a “boundary layer’’ effect in time or, more appropriately, an 
“initial layer’. The sense of approximation of the quasi-static theory and its 
relation to the dynamic theory have not, to the author’s knowledge, been previ- 
ously studied. 

In this paper we investigate these questions for a specific problem. It is the 
settling of a layer of a Kelvin material under a gravitational force. This is 
treated as a problem in the three-dimensional linear theory of viscoelasticity. 
We show, using a boundary layer expansion method [4]*, the manner in which 
the quasi-static solution approximates the dynamic solution for this problem. 
In addition, we obtain a sequence of approximating problems whose solutions 
give higher order or “‘dynamic’’ corrections to the quasi-static solution. Since 
the dynamic problem is easily solved, we can evaluate the accuracy of the 
approximations given by our expansion procedure. 


2. Formulation 


Consider an infinite layer of height y=H of a Voigt or Kelvin material 
resting on the x, z plane and bounded by smooth vertical walls at *»=-+L. 
At time ¢=0 we assume that gravity is ‘‘switched on’. We investigate the settling 
of this layer within the scope of the three-dimensional linear theory of dynamic 
viscoelasticity [6]. 

If o;;(%p, 4), 1,7, R=41, 2, 3**, are the components of stress and u,(x,, ¢) are 
the components of displacement, then the equations of motion are 


(1a) yy ge +F=o ous i= 2, 3. 


x Se ad s W i i 

; See also a related study by Wuirnam [5] concerning problems in wave propa- 
gation 
oP 


x* Here; 44=4, 44, Ves 
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Here g is the density of the material, and F,(x,, 2) are the components of the 
body force given by 


(1b) R90, b= ex 
where g is the gravitational constant. 
Stress-strain relations describing a Kelvin material are [6] 


, O65; 
en 


(2) 0;;=460;;+ 2ue,;; | ji ss 0; j | 2u 


where 0;; is Kronecker’s delta and e; ;(x,, ¢) are the components of strain. They 
are given by 


4 [ OU; Ou; 
a ée.= a 
Ba) 2 es = 
and 
Ou; 
b CS th 
(3) 2s 


The material constants A and w are the Lamé elastic constants, and 4’ and.’ 
are the corresponding viscous constants. It is frequently assumed that the 
Kelvin material is perfectly elastic with respect to dilatational motions. Then 
there are only three independent material constants. This assumption is not 
made here since the more general form (2) produces no additional complications. 


To complete the formulation appropriate boundary and initial conditions 
are required. We employ the following conditions: The surface y=H is force 
free, and therefore 


(4a) Gop 05 (Xda 2, t= O45 (Xd, 2b) =O 


The walls at x= +L are smooth and rigid so that 


(4b) 6,0(--L, y, 2, t) =0,5(4E, y, 2, t) =u (£L,y, Zt) =0. 


The material adheres to the x, z plane; hence 
(4c) WKF One id) = 05 t= Ag 253% 
Appropriate initial conditions are 


” _ Ou; (%, Y, 2,0) ae 
(5) u;(x, y, 2,0) =—— a O; a te 


This implies that the body is at rest at #=0*. The complete formulation of the 
dynamic viscoelasticity problem consists in equations (1—5). 


3. The Dynamic Problem 


We first introduce the dimensionless variables 


) Dyfi 
a ee a ZN ES SS (esarnke\ 
(6a) Uae seo ES Se Cacae 


* From (2), (3) and (5) it follows that the body is also stress free at /=0. 
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and the dimensionless material parameters 


2 gH*(A+2u) =, 74 (4 -). 
(6b) a Cera => O, a je yy 7 bu 
Then we seek a solution of (1—5) of the form 
yf ef Fe : 
(7) ty (%, ¥, 2, t) = Ug (%, 9,2, =0, Up(%, 4, 2,1) = Geer a (Vi Tes) 
We find that this is a solution if v(Y, T; «) is a solution of 
2 u au Sees 
(8a) av? + eyzar ~ © are?) 
at, dv (4, T; €) OFo (Aealege)\ ae 
(8b) v(0, T; e) =0, ae | ae 
GUY 0? 
(8c) OY, Oca) = é) = OP 


We further define dimensionless stress components as 


A+2 A+ 2 
o(Y, Tie) = (TERE Joule 28),  0,(¥, Ts) =( 7h osa(a, 9,20), 


oy (YN; sey é) cae (og H)* On (%, y; 4; t) $ 
Then (2), (6) and (7) yield 


a Ou mor av 
(9) %= a= sy tUtdaaer W=or+ ayer: 


All other stress components are identically zero. 
A formal solution of (8) is 


(10a) OS pres ig a Sry ey eer sin B,, Y. 


Here 67 — ee and 


(10b) v'(Y) =1(¥2—2Y) 


is the solution of the corresponding static elastic problem. If ¢ is restricted to 
the range 0<e<a/4, then 


(10c) eae. b, exp [a, T]—a, exp [bn T) . 
where by—ay, 
ca 2 4 
(10d) ea z i é : z] {1 1 (zee) | fs 
= 2 
eet area se 
Ps Fa OL 
"78 (2n—1)8 


are the Fourier sine coefficients of — v*(Y). We observe that for e in this range 
a, and b, are negative for all ». Thus as T increases the solution (10) decays 
to v°(Y). If e>a/4, then some /,(7) have an oscillating decay in T. From (9) 
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and (10) the dynamic viscoelastic stresses are 


oy(¥,T32)=o-+ >) AnBy [fa (T) + 207 | cos, ¥, 
(11a) nt a 

o,(¥,T; 6) =0,(¥, Tse) =o+ Y AnBu|fa(T) + (1+a) 7) Joos, y, 
where ae 
(11 b) o°(Y) =0%(Y) =o(Y)=Y—1 


are the dimensionless static elastic stresses. 


It can be shown that the infinite series in (10) and (11) converge uniformly 
and absolutely for OX Y<1 and T=0 and that term-by-term differentiation 


_ of (10a) is permissible for the orders indicated in (8a) for all T=6 where 6>0. 


Thus the formal solution (10) is a solution of (8). Uniqueness of the solution of 
(8) can be shown by an “energy” argument*. In Section 6 the solution given by 
(10), which we call the dynamic solution, is compared with the solutions de- 
termined by the quasi-static and higher order theories. 


4. Derivation of the Quasi-Static Problems 


The differential equation of the quasi-static theory is obtained from the 
dynamic theory (8a) by setting e=0. The resulting equation is of first order 
in time, and hence its solution cannot satisfy both initial conditions (8c). Thus 
there is a “‘boundary layer’’, or, more properly, an “initial layer’? behavior in 
the variable T near T=0. To investigate systematically this behavior and to 
determine the appropriate initial condition for the quasi-static equation we 
assume that v(Y, 7; «) can be represented asymptotically by a power series 
ie 6": 

(12) w(Y,) Dee) ~~) 0" (Ye Pee 

n=0 
Here the v"(Y, 7) are called the quasi-static displacement coefficients, and we 
define v'=0 if ~<0. Similar expansions are also assumed for each stress com- 
ponent in (9). Substituting these expressions into (8a, b) and (9) and equating 
coefficients of like powers of k yields the sequence of equations: 


a? u" GP ots ee 
(13 a) ey? ey? eT onl Gare 
a av", T) , 20% (4, T) _ 
3.2) ioe) oY aveT 9 
Ov o2u" 


o" (Y, T) =o" (Y,T) = 


2 


(13c) 
- “~ ov" o2u" 
oy (VY, I) = 35 + gyar’ 


* The author is indebted to Professor K. O. FRiEpDRIcHs for bringing this to his 
attention. 
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Each of the differential equations in (13a) is of first order in T and therefore 
requires only one initial condition. Thus the expansion (12) cannot, in general, 
satisfy both initial conditions (8c). If the series (12) represents the solution of 
the dynamic problem, it can do so only in a T-interval away from T=0. The 
expansion does not hold uniformly up to T=O and is not valid there. The 
solution of the dynamic theory deviates from (12) as T—0 (see also the results 
presented in Section 6). We call the region of sharp deviation the “initial layer’. 
In the limit as ¢—>0 the initial layer vanishes (see Section 6). To obtain an 
expansion of the solution near T=0 we introduce a new independent variable, 
t, by “‘stretching’’ T [4]. We do this so that as e->0 the initial layer in the tr 
variable does not vanish. Then the resulting expansions in t may uniformly 
represent the solution up to T=0. Let 


(14) te tie 


Then any fixed neighborhood of zero in the T variable is transformed by (14) 
into an arbitrarily large neighborhood in the 7 variable provided « is sufficiently 
small. 


We define the initial layer quantities, 


VYe es) OMY eee SS (Motte) ail Vou deel 
(15) 
LAY 302 8) =O, sand Ola | ee ree 


It is assumed that each of them can be represented by a power series in k= é?: 
(16) VAY ev elime Vel eves 
n=0 


where V"(Y,7)=0 if <0. Substituting (14—16) and the corresponding ex- 
pansions for S and & into (8) and (9) and equating coefficients of like powers 
of k, we obtain the following initial layer problems: 


ayn gyn. geyn-t 
he Ce ee oe. 
yoo, AU» OPH) Lg 

; Mbp | oY : 
(17) 

a a ar"(Y, 0) 
We y 30) seen) ew ee = 

(Y, 0) rE 0, 

at oVv"-1 e2 V” » ov n—1 o2 yn 
Sn | ii! (rs NA oe | to PORES ae 

ay oes ae ay tT Utes ys 70.1 ae 


In addition, we require that the initial layer coefficients approach or ‘match’ 
the quasi-static coefficients for large t. More specifically we assume that the 


v"(Y, 7) are regular functions of T near T=O and expand them in a Taylor 
series, 


yt (es T) » pn (eis) ee 
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where 
1 6™v"(Y, 0) Ss 
(18a) pia lml ere 7 n=” 
OF are 

We find from (12) and (14) that near T=0 
(18b) v(¥, Tie) ~Y Dos(vy Te => UY, aR, 

n=0 m=0 n=0 
where 
(18¢) CYT) Dep (Y ee 

m=0 


are the quasi-static displacement coefficients near T=0 as functions of Y and t. 
Employing (15), (16) and (18), the matching conditions or the asymptotic forms 
for the V” are obtained by identifying each V” with the corresponding U” as 
é—>0. Using (14), we write these conditions as 

(18d) lime [V" (Yea) (YY, 7) = 0, T= 0 Aare 


TOC 


‘For example, by first setting »=0 and then 1 in (18) we find that 


Lo) Sin, (yO im 


tT >©CO TOC Te 


Equations (19) yield the appropriate initial conditions for v°(Y, 7) and v!(Y, 7) 
in terms of the asymptotic form of the solutions of the initial layer problems, 
(17), of order zero and one. Appropriate initial conditions for v"(Y, 7) with 
n> are obtained in a similar manner from (18). 


5. The Asymptotic Form of V(Y, Tt) 


To obtain, from (18), the proper initial values for the quasi-static problems 
(13) we require the asymptotic forms in t of V”. These are obtained from the 
solutions of the initial layer problems (17). The zero™ order initial layer problem 
is given by (17) with* 7=0 and has the solution 


(20) VU Nt) Oe 
Therefore, from the first of (19) the appropriate initial condition for v°(Y, T) is 
(21) v°(Y, 0) =O. 

The first-order initial layer problem is given by (20) and (17) with »=1: 


a V1 Gaee 


oY20t .+=«OT? : 
é2V1(4, 7) 
if Ee zat 
(22) V*(0, 1) oY at 0, 
éV1(Y, 0) 
71 — = — 
V1(Y, 0) 0. 


* We recall that V"=R"=S"=0 if n<0. 
Arch. Rational Mech. Anal., Vol. 7 27 
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By employing Laplace transform techniques, we find that an integral represen- 
tation of the solution of (22) is 


esa review petras, 
where , 
(23 b) G (Veg a ene ae 


p* cosh p® 


and c is chosen so that the path of integration is to the right of the singularities 
of G(Y, p). The integral in (23) can be evaluated in terms of theta functions 
[7], yielding an approximation to the dynamic solution near T=0. However, 
to obtain from the second of (19) the correct initial condition for v1(Y, T) we 
require V1(Y, 1) for large t. This is easily evaluated from (23) by standard 
methods [8]. Substituting the result of this calculation into the second of (19), 
we obtain the initial condition for v1(Y, 7) as 


(24) v(Y,0) =2, [¥*— 4¥? + 8Y], 
and the auxiliary condition, 
CUTTY, © far ae : 
ee (Y2— 2Y). 


The latter condition is automatically satisfied; see (25a) below. 


The initial conditions (21) and (24) complete the formulation of the quasi- 
static problems of the two lowest orders. The “‘zero™ order quasi-static problem’, 
1.e. the quasi-static theory consisting of (21) and (13) with »=0, has the solution 


(25a) v(Y, T) =4(¥?—2Y)(1—e7-7), 
(28D) oY lye (Y = 4) ao CY hy =a y 2) ee ee eae 


The “first-order quasi-static problem” consisting in (24) and (13) with n=1 
has the solution 


(26a) WY, T)=—f(Y) = Te-?, 
(26b) of(¥, T)=BM 7, Al, T)=oh(¥, T) =o (Y, Tt —a(T —2)], 
where 
VAN eye = 
(26c) (Y=—4 pice Heel 


6. Discussion of Results 


The first approximation, valid away from T= 0, to the solution of the dynamic 
problem, (10) and (11), is given by (25). We observe from (9), (40b) and (41 b) 
that for the static elastic problem 

O11 


—_/ ¢ sas 
eEy: 22 22+ 
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The viscoelastic stress o$, equation (25b), “jumps” at T=0 to the static elastic 
value. Thus, within the first approximation o, shows no viscoelastic effect. 
The behavior of of and of depend strongly upon «; see Fig. 1. For <=0*, o8=o!, 
and there are no viscoelastic effects. However, for «>0, 0? “jumps” at T=0 
to of=(Y—1)(1+«) and then approaches its static elastic value from above: 


Dynamic Sol o,{0,7,€) 


— — — asi- Static Sol of(0,7) 


ares 


Fig. 1. The T-variation for fixed e and varying « of o, obtained from the dynamic and quasi-static solutions 


for «<0 the approach is from below. Therefore, according to the quasi-static 
theory for «>0, o2>0% and of is the maximum stress. This is in contrast to 
the static elastic case. In Fig. 1, 62 (0, 7) and the corresponding dynamic stress 
are given as functions of T for e=a/5 and varying «. From this figure we observe 
that for «> 0 the accuracy of approximation of o2 (0, T) decreases as T decreases. 
When T< 1.5, o° is, in general, a poor approximation. For e=2/10 (not shown) 
the quasi-static solution is an accurate approximation for (> .25 and —1<¢<1. 


* Te. N= 4. 
le 


bo 
NI 
* 
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For ¢=2/100 (not shown) the approximation is accurate for T> O01 and the 
same range of «. In Figs. 2 and 3 the quasi-static and dynamic solutions are 
compared for fixed « and varying «. The vanishing of the initial layer as e—>0 
is illustrated in these figures. 


20 


Dynamic Sol. 0,(0,7;€) 
76 


— — — uasi- Static Sol 6f(0,7) 


G27 


76 


74 


72 


707} 


- 0; (0, T,€) 


08 


06 


02 


i 13 | | 
% 7 a: Bi > 5 i 6 


Fig. 2. The T-variation for fixed « and varying ¢ of oz obtained from the dynamic and quasi-static solutions 


To obtain more accurate approximations for “larger”? ¢« and T away from 


zero, the solutions of the zero‘ and first order quasi-static problems, (24, 25), 
are combined according to (12) to form: 


(iy a 0 eee | QQ) eae: 2 
Oy =G,+&°O,, Go, =G6, =, + &70,. 


In Fig. 3, the improved accuracy of of? as compared with o° is apparent. 


More accurate approximations to the dynamic solution for T near zero may 
be obtained from the solutions of the initial layer problems, (17). 


In conclusion we note that our approximation procedure can also be obtained 
directly from the dynamic solution. Our method is advantageous for those 
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problems whose dynamic solution cannot be found explicitly. Increasingly 
accurate approximations may then be obtained* by solving a sequence of com- 
paratively simple problems. 


e=70 ,andoy (0,7) ond oj (0,7; = 


1 100 


-0, (0,7, €) 


a=0 
Dynamic Sol. 00,7; é) 


SS ofa,fe) 


% 7 2 of 4 a 7 6 


Fig. 3. The 7-variation of oy, obtained from the dynamic, the quasi-static and the “improved” quasi-static solution 


The work reported in this paper was performed under Contract No. AT(30-1)-1480 
with the United States Atomic Energy Commission. The author is indebted to 
H. B. KeLver for his many helpful suggestions and stimulating discussions. 


* Provided ¢ is sufficiently small. 
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For ¢=2/100 (not shown) the approximation is accurate for T> .01 and the 
same range of «. In Figs. 2 and 3 the quasi-static and dynamic solutions are 
compared for fixed « and varying «. The vanishing of the initial layer as e>0 
is illustrated in these figures. 


Dynamite Sol. 0,(0,7;€) 


— — — uasi-Static Sol. of (0,7) 


Ort 


| 
“9 7 “a 3 4 on at ne 


Fig. 2. The T-variation for fixed « and varying ¢ of oz obtained from the dynamic and quasi-static solutions 


To obtain more accurate approximations for “‘larger’’ ¢« and T away from 
zero, the solutions of the zero'* and first order quasi-static problems, (24, 25), 
are combined according to (12) to form: 


oP = 08+ stot, o® =o! <0 + etal. 
In Fig. 3, the improved accuracy of of? as compared with o}, is apparent. 
More accurate approximations to the dynamic solution for T near zero may 
be obtained from the solutions of the initial layer problems, (17). 


In conclusion we note that our approximation procedure can also be obtained 
directly from the dynamic solution. Our method is advantageous for those 
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accurate approximations may then be obtained* by solving a sequence of com- 
paratively simple problems. 


700 


‘€= ,and oy (0,7) and a (0,7; 
4 


- 0, (0,7; €) 


a=0 
Dynamic Sol. oy (0,7;€) 


Se Oe, 


% 7 2 ow ¥ oF 7 6 


Fig. 3. The T-variation of oy obtained from the dynamic, the quasi-static and the “improved” quasi-static solution 
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with the United States Atomic Energy Commission. The author is indebted to 
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Eine Verscharfung der asymptotischen Gesetze 


eleRtromagnetischer Hoblraumschwingungen 


Horst NIEMEYER 


Vorgelegt von C. MULLER 


Einleitung 


Es sei G ein regulares! Gebiet des dreidimensionalen euklidischen Raumes, 
das mit homogener und isotroper Materie erfiillt ist und das von einer im physi- 
kalischen Sinne ideal leitenden, dreimal stetig differenzierbaren Randflache F 
begrenzt wird. Die Materialeigenschaften seien durch die Dielektrizitatskon- 
stante « und die Permeabilitat ~ festgelegt. Die Konstanten e und yw seien 
positiv reell. 


Elektromagnetische Eigenschwingungen ©, $ zur Eigenfrequenz w geniigen 
in G den Gleichungen 


Vx$+¢mweE=0, VxE—couhH=o 


und der Randbedingung nx€=p auf F. Dabei ist n der ins AuBere von G 
weisende Normalenvektor auf &. Mathematisch stellt dies ein Eigenwertproblem 
fiir ein System von sechs partiellen Differentialgleichungen fiir die insgesamt 
sechs Komponentenfunktionen der Vektoren © und § dar mit einer homogenen 
Randbedingung fiir die Komponenten von 6. 


In einer friiheren Arbeit gemeinsam mit C. MULLER [9] wurde gezeigt, daB 
dieses Eigenwertproblem unendlich viele Eigenfrequenzen @,, @:, @3, ... besitzt, 
die alle reell sind, sich im Endlichen nicht haufen und jeweils von endlicher 
Vielfachheit sind. Der Fall w=0 ist hierbei auszuschlieBen. Die zugehérigen 
Figenschwingungen ©, und §,, sind in G stetig differenzierbar, im abgeschlossenen 
Bereich G stetig und erfiillen die Randbedingung nx ©,=o auf F. Jeder Eigen- 
frequenz , entspricht eine Eigenfrequenz —w,, gleicher Vielfachheit, wobei die 
Eigenschwingungen zur Eigenfrequenz —q,, konjugiert komplex zu denen von 
@, sind. Wir betrachten daher nur die positiven Eigenfrequenzen, die wir uns 
monoton wachsend angeordnet denken wollen (0<@,S@,S:::). 


In der erwahnten Arbeit wurden fiir die so normierten Eigenschwingungen 
und Eigenfrequenzen unter anderem folgende asymptotische Aussagen bewiesen 


? , regular im Sinne von O. D. KEttoc [6], S. 113. 
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CEG, SOWs « .\/2 
("EG, x00) > IG)? =, (ew)? *, + 0(x) 


eRe wr 
- 3 
» [6 @l? =n (ew? =, + o (x9) 
0<anS% we 
V- x3 
Nec! EAD ae ct) 
<OnS% 


Hierbei ist V das Volumen des betrachteten Gebietes G, und >» 1 bedeutet 


° fas . ° . 9<onS* 
die Anzahl der positiven Eigenfrequenzen ,, die kleiner oder gleich x sind. 


Wir werden in dieser Arbeit zeigen, daB die Fehlerglieder in diesen asym- 
ptotischen Aussagen verbessert werden kénnen, und zwar beweisen wir (r€G, 


) €G) 
©, (x)|? 3 
yl ee =p(e w)t= + 0(1) 
0<a,S% _ 
vas 
>» 16 @)? =HleniA, +0(x) 
0<anS% sn 
0<ayn,Sx se 
und 
1 = (eu)? or + O (x? log x). 
0<anS%x 


In der vorletzten Gleichung ist dabei der als Spalte geschriebene Vektor 6, (r) 
mit dem transponierten Vektor E>(y) nach den Regeln des Matrizenkalkiils zu 
multiplizierent. In der entstehenden Matrizensumme bleibt dann jede Kompo- 
nente beschrankt fiir 7—> oo. 

Asymptotische Aussagen iiber die Eigenwerte der Helmholtzschen Schwin- 
gungsgleichung sind zuerst von H. WeytL bewiesen worden [10], [J/], [12]. Im 
Jahre 1934 konnte T. CARLEMAN [3], [4] dariiber hinaus auch fiir die Eigen- 
funktionen asymptotische Beziehungen herleiten. 

Sind uw, die zur Randbedingung w=0 auf F gehérigen orthonormierten Eigen- 
funktionen der (dreidimensionalen) Schwingungsgleichung und 0<A,S/,S:-- 
die zugehérigen Eigenwerte, so gilt nach T. CARLEMAN 


Un Ve i 
7 Re py roe te oN 
n=* 3 
»y “w= me + 0 (x?) 
An ies 

SS V x! 5 

> tages 672 stags) 
Anat 


Die letztgenannte Abschatzung der Anzahl der Eigenwerte ist die ursprting- 
liche Weylsche Formel, die spater von R. Courant mit Methoden der Variations- 
rechnung verbessert wurde. CouRANT konnte zeigen, da das Fehlerglied in 
dieser Formel nur von der Ordnung O(x - log x) ist. 


1 Transponieren deuten wir durch ~ an, den Ubergang zum konjugiert komplexen 
durch ~. 
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V. G. Avakumovié [1], [2] verscharfte 1952 auch die asymptotischen Be- 
ziehungen iiber die Eigenfunktionen und bewies: 


> ORE, Se VF OVA} 


V Digi 
AnS*% 
y Unt he () 04) (e+) 
InS : 
y ae = jax +O(2). 


Entsprechende Relationen sind von AvAkuMovic¢ auch fiir die Eigenfunktionen 
der Schwingungsgleichung auf geschlossenen Riemannschen Mannigfaltigkeiten 
abgeleitet worden. Im Fall geschlossener Riemannscher Mannigfaltigkeiten gelang 
es ihm auch, durch Integration der letzten ree fiir die Anzahl der Eigen- 
werte die Formel 

3 ‘— Ae O (x) 

AInS% 
zu erhalten. Im Falle berandeter Mannigfaltigkeiten ist jedoch ein besseres als 
das Courantsche Ergebnis bis jetzt nicht bekannt. 

Wir zeigen hier, daB entsprechende verscharfte Aussagen fiir die elektro- 
magnetischen Eigenschwingungen gelten. Der Beweisgedanke, der in ahnlicher 
Weise von T. CARLEMAN und V. G. AVAKUMOVIC¢ angewandt worden ist, geht von 
Bilinearentwicklungen aus, die die Eigenschwingungen mit den Greenschen Funk- 
tionen bzw. Greenschen Tensoren verkniipfen. Im Fall der Schwingungsgleichung 
kann man z.B. von gtr sniveaaai 


ae Sr ee ) = G(r, 930) — G9; — 7) 


ausgehen, wobei G(r, ); — t?) die Greensche Funktion der Schwingungsgleichung 
(A —t?)u=0 mit der Randbedingung w=0 auf F ist. 


Setzen wir HE SINS 
O(x) = <= ? 
ea An 
so ist () oo 
Un(£) Un(Q) os dO(x) _ Ren 2 


0 
Wir kénnen da her die Bilinearentwicklung als Stieltjes-Transformation der Funk- 
tion @(x) auffassen. Wollen wir aus dem Verhalten der Greenschen Funktion 
fiir t->oo auf das Verhalten von O(x) fiir x->co schlieBen, so fiihrt das auf einen 
Fragenkreis, der durch den Begriff ,,Taubersche Satze‘‘ gekennzeichnet ist. Wir 
werden hier einen Tauberschen Satz fiir die Laplace-Transformationen benutzen, 
der von A. E. IncHaAm [5] 1935 bewiesen worden ist. 


In §1 stellen wir die wichtigsten Eigenschaften der Greenschen Tensoren fiir 
die Maxwellschen Gleichungen zusammen, die wir spaterhin benétigen werden. 
In §2 wird das Verhalten der Bilinearentwicklung bzw. der Greenschen Tensoren 


fiir t—>oo untersucht und in §§3 und 4 wird der erwahnte Taubersche Satz auf 
die Bilinearentwicklung angewandt. 
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§ 1. Greensche Tensoren 


In diesem Paragraphen stellen wir die wichtigsten Eigenschaften der Green- 
schen Tensoren fiir die Maxwellschen Gleichungen zusammen, die wir spater 
benétigen. Die Ergebnisse sind in der erwahnten Arbeit [9] enthalten. 

Es sei G ein endliches, reguléres Gebiet, das von der orientierbaren, zusammen- 
hangenden, dreimal stetig differenzierbaren Randfliche F berandet wird. In G 
betrachten wir das System der Maxwellschen Gleichung 


(14) VxD+¢weE=0, VxE—¢wouH=o 


mit der Randbedingung nx €=p auf F. Wir setzen ¢ und y als positiv reell 
voraus. 

Wie bereits bemerkt, gibt es unendlich viele (reelle) Eigenfrequenzen w,, und 
zugehorige Eigenschwingungen ©,, $,, wobei wir uns auf die positiven Eigen- 
frequenzen beschranken und diese monoton wachsend anordnen wollen. Sind ¢ 
und mw von 1 verschieden, so kénnen wir die Gleichungen (1.1) durch die Sub- 


stitutionen 


(1.2) H=yeh, C=eE, ofen=a 
Zu 
(1.3) VxE—¢oH=0, Vxh+¢o€=0 


mit der Randbedingung n x =p auf F normieren. Wir setzen daher weiterhin 
€=u=1 voraus. 

Greensche Tensoren &,;, (x, ); @) (¢, R=1, 2) zu den Maxwellschen Gleichungen 
k6énnen durch folgende Relationen definiert werden: 


1. Ihre singularen Anteile &%,(r, ); ) sind durch die Forderung bestimmt, 


daB fiir jeden konstanten Vektor a und mit @(r, y)= — 
OT (z, 9; w)a=¢oPa+ <P, (a V, ) 
G(x, ); o)a=V,Oxa 

(1.4) 
Gis (x,y; @)a=—V,®xa 
Gia (t, 9; @)a=cw a+ V,(aV, ©) 

ist. 

2. Die kompensierenden Anteile 
(1.5) Rp (U0; ©) = Gia, 9; @) — Gee, 9; @) 


sind in jeder Komponente als Funktion von yx stetig differenzierbar in G. 


3. Die Felder G,,a und G,,a (k=1, 2) erfiillen als Funktionen von x fir 
y-+ty die homogenen Maxwellschen Gleichungen, d.h. 
V, XG 1,4 — ¢@ G,a=0 
ue) Vx G2,a+ ¢o O,,a=0. 
4. Es gelten die Randbedingungen 


(1.7) n(c')X@Gi,(e',9;@)a=o (EF, y€G). 
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Ist w keine Eigenfrequenz der Maxwellschen Gleichung — insbesondere also fiir 
rein imaginire Werte von m —, so sind die Greenschen Tensoren eindeutig 
bestimmt. Der Existenznachweis der Greenschen Tensoren fiihrt auf ein Rand- 
wertproblem der Maxwellschen Gleichungen, das von C. MULLER [7], [8] bereits 
gelost worden ist. 

Fiir die kompensierenden Anteile macht man hierzu mit einem noch unbe- 
kannten Tangentialfeld v(p, y) an F (p CF, )€G) den Ansatz 
(1.8) (2,9; 0)a=—{- f v(p,n) x% Olp,x) dF. 

F 


Die Randbedingung Gl. (1.7) liefert dann nach den Sprungrelationen die Integral- 
gleichung 


(1.9) (py) — 5 fm (e) x (0(4, 9) x7, (a, p)) €F = 2n(p) x [Mp Pp, 0) Xa] 


symbolisch 
(1.10) b—Mv=w(p,y) (= 2n(p) x VY, ® xa)) 
geschrieben. 

Fiihren wir als Norm des Flachenfeldes » (p, \) 


|» | = supr |v (p, »)| 
DEY 


ein, so ist die Integraltransformation M vollstetig im Banach-Raum aller stetigen 
Flachenfelder auf F’. Ist w keine Eigenfrequenz der Gl. (3), so ist die Integral- 
gleichung (10) fiir beliebig stetige Flachenfelder eindeutig lésbar. Dieselbe Inte- 
gralgleichung mit anderem inhomogenen stetigen Flachenfeld erhalt man fiir das 
erzeugende Feld des kompensierenden Anteils &,, a. 

Die mit den Losungen dieser Integralgleichung gewonnenen Tensoren 
M1 (z, ); @) und R(x, y; @), sowie die daraus durch Differentiation gewonnenen 
Tensoren Sy; (r, 9; @) und &,,(z, ); w) erfiillen alle Forderungen in Gl. (1.4) bis 
(Gear 

Fiir die Greenschen Tensoren gelten — ebenso wie fiir die singularen und 
kompensierenden Anteile allein — folgende Symmetrierelationen 


(1.11) On (zy, 9; @) = ©, (b, 23 w), 
(4:42) Gro (vt, 9; wo) = — Ge, (y, x; @), 
(1.13) 22 (x, 0; w) = Goo(y, 2; @). 


SchheBlich kann man zeigen, daB die Tensoren ®yo(x%, 9; @) und G(r, y; w) 
auch noch — in einem gewissen ausgearteten Sinn — fiir w=0 existieren. 


Bildet man nun mit den Tensoren ®,, und @,, folgende iterierte Tensoren 


(1.14) Gz (vt, 9; w) = — —s i. yo (v, 3; 0) G1 (3, 9; @) adV;, 
G 


(1.15) Gat, 93 @) = — Teg Sor (ts §5 0) G12, 95 @) ah, 
G 
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so erfiillen die elektromagnetischen Eigenschwingungen die Integralgleichungen 


(1.16) (wo, + 7?) f Ge(e, 9; 27) E, (y) dV, = G, (x), 
G 

(1.17) (ow, + 2) f Gy (z, 9; £7) 8, (b) dV = 9, (2), 
G 


und es gelten die Bilinearentwicklungen 


Mees . En (x ) GK (y) 
(1.18) Gr (x, 9; 0) — Gee, 9; £7) a aay oh (ozt 18)? 

ie he Hn (z) ie 
(1.19) Gx, 9; 0) — Gy(e, 9; 47) a “wo? (wo? 


§ 2. Die Bilinearentwicklung 


Den Ausgangspunkt unserer Betrachtungen bildet die Bilinearentwicklung 
Gl. (1.18) 


§ 
(2.1) 2 Y Se BO) = Gp (¢,; 0) — Gale, 93 0). 


bas 0- 


Lassen wir )—>r streben und bilden die Spur der Matrizen auf beiden Seiten, 
so ergibt sich 


(2.2) Lead ote wey me Spur {Gz (x, ); 0) — Ge (x, 9; ¢7)}. 


osot, 


Wir wollen das Verhalten der rechten Seite fiir t—+oco untersuchen. Dazu zer- 
legen wir die Tensoren nach Gl. (1.4) und (1.5) in ihre singularen und kompen- 
sierenden Bestandteile und erhalten drei Summanden, die wir mit R,, R, und Rk, 
bezeichnen wollen: 


(2.3) Spur {Gz (z, 9; 0) — Gelv, 93 ¢7)} = R(t) + Rye, t) + R(t), 


yr 
mit 
CO es(An) hae) li poy i) Gts (x, 3; 0) {G21 (3, 3 47) — G3, 9; O)$ AV, 
>t 
(421)? Re (x, t) = Spur f G2 (x, 3; 0) Ky (5,45 27) A+ 
(2.5) Ge : 
+ Spur f &12(t, 5; 0) Gai (g, 83 47) aVy 
G 
und 
— (42)? g(t) = Spur f &y2(t, 35 0) Goi (g, 25 0) dV, + 
(2.6) oe ; ; 
+ Spur f Gf2(x, 3; 0) Kei (5,25 0) dV. 
G 
Wir bemerken, daB R(x) — wie schon in der Schreibweise angedeutet — 


nicht von t abhangt. Daher diskutieren wir zunachst das Verhalten von hk, 
und Ry. 
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Lemma 1. Fiiy r€G, v,= inf |p—z|2a>0! und t=O ist 
os 


(2.7) RG) =s-+Ge) +0(Fe*)). 
Die Funktion C,(r) hangt dabei nicht von +t ab und geniigt der Abschatzung 
|. @)| = 2 (A/a) 

Beweis. Nach Gl. (4.4) ist 


(2.8) Spur Gre (x, 3; 0) GF (3, 9; 47) = — 2V, a -V; D3, 9), 
und daher 
(2.9) eg s Spur f Of 1,3; 0) OF, (3, 9; a ax | Bia (426.0) 4%, 


Den letzten Ausdruck formen wir mit Hilfe des GauBschen Integralsatzes um 
und erhalten 


10) [Pama V, Bis, 9) dV, = 42 Ee, 9) + f OG, 0) aR, 


éng |3—2| 
F 
fig 
1 1 7 
(2.11) Ry (ot) = Se f( D3, »)) 5 
F 


ls—zl és |3—tl 
Fir alle Punkte r€G mit 7,2«>0 und fir alle r20 ist 
| ae] 1 eWet id 1 her ee 
Dy é a1 bdo Ovens [| 2 dF =O|= ). 
ee if [t=a), Ore laa) ( x J om least * ise 
F F 
Ferner ist 


1 a) 1 
os : 4 = Git 
a) I tea eng le-al 7 1(¢) 


eine von t unabhangige Funktion von y, fiir die entsprechend die Abschatzung 
|C, (x)| =O (1/«) gilt. Damit ist Lemma 14 bewiesen. 


Wir zeigen weiterhin 
Lemma 2. Fiir r€G, 1,=a>0 und tZ0 18¢ 
a/2 2 
(2.14) R,(t, t) = few. rye" (14+cr)dr +0(4 e783 Me 
0 
Die Funktion g(r, x) ist dabei unabhdngig von t, bei festem x als Funktion von r 
im Intervall 07S «a/2 stetig, und geniigt dort gleichmafig der Abschitzung O(1 hed 


Beweis. a) Nach Gl. (1.8) und der Symmetrierelation Gl. (1.12) erhalten wir 
den Tensor S.; (3,4; 47) mit Hilfe eines Flachenfeldes v(p, 3) (p CF) gemaB2 


(2.15) PE ts A) at sa aN aba Te P(p, x) dF 


1 Es ist also 7, der Abstand des Punktes r von der Randflache F. 


® Wir wollen mit x, ), 3 Punkte aus G, mit p, q Punkte auf F und mit pv, ty Tangential- 
felder an F bezeichnen. 
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wobei » die eindeutig bestimmte Lésung der Integralgleichung 


(2.16) b — Mv =w/(p, 3) = 2n(p) x [V, Dp, 3) xa] 
mit 

(2.17) yi errraifil has » (a, 8) x V7q P(a, b)) aE 
ist. 


b) Das Flachenfeld tw (p, 3) wird singular fiir 3->p, und zwar gilt 


ie = 1 
(2.18) [ (, 8)| = 0 (| P(p,a)|) =O( 5) 
und 
(2.19) =O(1) (unabhangig von @ und 4). 
is—plSe 


Beide Abschatzungen gelten gleichmaBig beziiglich r=0. Wir kénnen daher ein 
Ergebnis einer friiheren Arbeit [9] verwenden (Lemma 6 und 8 dort), nach dem 


1 
Mw] = ] 
Biel Os 
2.20 M?tp| = O10 Leo 
oe) | | s rani 
[M3 | = 
ist, und zwar ebenfalls gleichmaBig beziiglich aller t= 0. 
c) Fiihren wir jetzt das Flachenfeld 
(2.21) b* =b — vw — Mw — M* 


ein, so erfiillt »* die Integralgleichung 


(2.22) b* — Mv* = Mp = w*. 


Nun konvergiert die Neumannsche Reihe dieser Integralgleichung? fiir geniigend 
groBe t, denn es ist mit 7=|p —q| 


Ions a J |r) x (0*(q. 8) xVa - | e V4 7) af, 
(2.23) =o feeruten af) [| 
F 
=|p*]-0(-). 
Es gibt daher eine nur vom Gebiet G abhangige Konstante 1,, so daf 
|, v*] S 2 [o*| 


und somit nach Gl. (22) 
(2.24) |o*| Ss 2 m*| 


ist fiir alle r= Tp. 


1 Wir bezeichnen die Abhangigkeit des Operators M von t durch /,. 
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d) Fiir alle tr mit OX tS, ist die Integralgleichung (24) eindeutig lésbar. 
Fiir t=0 ergibt sich dies aus [9], Lemma 5. Der Operator M, hangt auBerdem 
stetig von t ab, d.h. es ist 


(2.25) (4, — M,,) 0*| S D,| 1 — 7| - |v* 


, 


mit einer nur vom Gebiet G abhangigen Konstanten D,. Es gibt daher eine 
Konstante D,, so da fiir alle t mit O05 TST, 


(2.26) |v*] S D,- | w*| 
ist. Andernfalls namlich gabe es Folgen t,, (0< 7, ST), 


be (lor f) md wy 
mit lim ||w*]]=0 und 
n—>oo 


Th 


(2.27) De Goisaine: 


Wir kénnen (gegebenenfalls nach Bildung von Teilfolgen) annehmen, daB 
lim 7,=t* ist (0O=7t*=7)). Dann ware 


n—> CO 


(2.28) y* — M..o* = w* + (M, — M,) oF. 


Tn %% 


Da die Integralgleichung fiir r=t* eindeutig aufldsbar ist, gilt mit einer Kon- 
stanten D, 
lon | SDs |e + (M,, — Mz) vp 


(2.29) 4 
< D;||w*|| + Ds D,- |t,, — t*]. 
Es ware demnach im. |v; =0 im Widerspruch zur Annahme |\p*||=1. 
e) Setzt man D= Max (Dg, 2), so ist 
(2.30) Jo*(p. 3) 


fiir alle 3 aus dem abgeschlossenen Bereich G und fiir alle r=0. Es folgt also 
nach Gl. (20) und (21) 
[vo] S|] +] Mw] +] M2w] +011), 


<D-|w*|=D-|M3 wl =0(1) 


oder 
(2.31) |v (p,3)] =O(|p —3 


gleichmaBig fiir alle r=0. 


) 


f) Aus Gl. (15) und (314) ergibt sich als Abschatzung! fiir den Tensor Ra 


> po tipal . 
| a1 25 ¢) =o padres eae in) 
F 


~ [p—rl? [pal 


=o[e# fip—al*-Ip—altan] 


F 


(2.32) 


fiir 7,=a>0 und alle r=0. 


1 Unter |S, ,| wollen wir die Summe der Absolutbetrage der einzelnen Komponenten 
von §,, verstehen. 


p 
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Wir erhalten nunmehr fiir den ersten Summanden in Gl. (5) 


- Spur | is (x, 5; 0) Stan Ges 2 U) dV, 


ia ae : : Fav. 
( J fat) boar aré var) 
dee 


ae, 
(4 2)? 


g) Es bleibt noch der zweite Summand in Gl. (5) zu betrachten: 
1 

Te. ok: f S235 0) Goi (5,05 7) aV; 

(2.34) 


= Gar =a fl Tis Or dV; i fi Ryo Roy an : 
& G 


Fiir den Tensor (rx, 3; 0) gilt die Abschatzung 


(2.35) [fats O(n), 
wobei Ry, = inf (|p —z|+|P —al) der ,,Lichtweg‘’ zwischen y und 4 ist ([9] 
p a 


Satz 6). Erst recht ist daher |®,.(x, 3; 0)| ot und wir erhalten mithin 
fiir den zweiten Anteil auf der rechten Seite von Gl. (34) wie in Gl. (33) die 


Abschatzung o(s e 2’). Bei dem ersten Anteil zerlegen wir das Integrations- 
a 


gebiet in zwei Teile, namlich in die Kugel G.={s, lr—3|S = und das Rest- 
gebiet G,=G—G,. 4 


Im Restgebiet G, ist |r te) d.h. 
eae 
(2.36) | BH (3, =0(e 2 Sienna 
eck = Py (a) 
(2.37) 


Cie: 
In der Kugel G, ist unter Benutzung von Gl. (1.4) mit r=|t—4| 


as Spur Ry 2 (£, §; 0) OF (3,45 7) AV, 


(2.38) Cs Af (Set one sa = 


ones 


Dabei sind e¢,, ¢, und e, die Einheitsvektoren in Richtung der Koordinatenachsen. 
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Zerlegen wir die Volumenintegration in eine Flachenintegration und eine 
Integration tiber den Kugelradius, so folgt 


a/2 


(2.39) te Spur [%.6h4aK= | ee, rhe "(A+tr) dr 
; 0 
mit 
; 3 
(2.40) gnu = (ant ii px Ro 1 (3, £5 0) &; ei) (V,7) af, 
le—al=r 


Da 8, (3,2; 0) in G, stetig ist als Funktion von 3, so ist auch g(r, zx) stetig im 
Intervall 0<r<a/2 als Funktion von 7. Uberdies geniigt g(7, x) nach Gl. (2.35) 
der Abschatzung 


(2.44) rolachentag ae ee 


lr—3l=r & 


und zwar gleichmabig fiir alle y mit 0X7 a@/2. Aus Gl. (33), (37) und (39) ergibt 
sich die Behauptung von Lemma 2. 


Wir kénnen nunmehr die in Gl. (2.2) auftretende Reihe abschatzen. Es gilt 


Satz 1. Fir r€G, 7,2a>0 und t20 ist 


2 |G, (x) |? fi (Le =F 2") 
(2.42) pas ORE ice r+ | gl xe "(14+ drto ; 
Die Funktion C(x) und g(r, x) hangen nicht von t ab. C(x) geniigt der Abschdtzung 
O(1/a), g(7, x) tst stetig als Funktion von r im Intervall 0SrSa/2, und es gilt 
[g(, 2)| =O (1/0?). 

Beweis. Nach Gl. (2.3) ist die in Gl. (42) linker Hand auftretende Reihe 
gleich R,(z, t) + Ry (rv, t) + Rs (x). Fiir R, (x, t) und R,(r, t) haben wir bereits in 
Lemma 1 und 2 gezeigt, daB sie die Behauptungen unseres Satzes erfiillen. Wir 


haben daher lediglich zu zeigen, daB die von t nicht abhangige Funktion 
R3(t)=O(1/a) ist. Nach Gl. (2.6) und (2.35) ist 


(2.43) 10 fae ea an), 
v3 


wobei Re) inf (fp +|p—s|) der Lichtweg zwischen x und 3 ist. Zerlegen 


wir das Intepradionseetiet wie im Beweis zu Lemma 2 (g) in die Kugel 
G,= {3, | —3| Sa/2} und das Restgebiet G.=G — G,, so ist in G,: R,, 2 $0, d-h. 


(2.44) 
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Im Restgebiet G, ist |r —3|=«/2, also 


2.4 freer th [pr en=o(4). 
ee) py ey NE Atl coer tu dl Een Sh ge a. 
Insgesamt ist also 

(2.46) Rel =0(2), 


und Satz 1 ist damit bewiesen. 


Zum AbschluB dieses Paragraphen beweisen wir eine zu Satz 1 analoge Ab- 
schatzung der Bilinearentwicklung Gl. (2.1) fiir y=—y. Wir nehmen dabei an, 
daB r€G, »E€G und y+ y ist, und setzen 


ty=a>0, |¢—y|=8>0, Min(@,p)=y>0. 


Wir bezeichnen die Komponenten des Vektors 6,,(r) mit E)(r), E% (x), E3 (zx). 
Dann gilt 


Satz 2. Fir y>0 und tZ0 1st 


ik (+) Fl ue ears 
(2.47) 7 >! ae = Cait, n+ f eeilrse, ye (+ t7) dr + Ole : ik 
On OMe ie 


0 


Dabet hiingen Cy, (x, 9) und gy(7; x, ) nicht von t ab. Die Funktionen g,,(r; t, Y) 
sind stetig im Intervall OSrSy/2 (k, 1=1, 2, 3). 


Beweis. Es ist nach Gl. (2.1) und Gl. (1.14) 


(2.48) Ie £0) Ent). — oF (Wp (2, 9; 0) — Oe, 95 2) ep 
duh. = 
k Tl 
ot) 2p) En) — — A f oF Gyalts95 0) Gaile 95 0) dK + 
On>0 a " 
(2.49) 


G 
+ aap | 8 » Gyo (t, 33 0) Gor (3, 9; 27) e: aVj. 


Nun setzen wir 


(2.50) Ceilt, 9) = — age f ee Gyo (L, 3; 0) Gar, 9; 0) ep dV. 
G 
Ferner ist nach Gl. (32) 


amr |  Sralt a5 0) or 9; 42) eat, —o(e * f lp a] *logy;*ar) 
‘ G 


(2.54) i =0(e#") =0(e"#"), 


In dem verbleibenden Anteil 


aaa | Bro (t, 55 0) O21 (3, 9; #7) dV; 
G 
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zerlegen wir das Integrationsgebiet wieder in die Kugel G,={3, |» —3|Sy/2} 
und das Restgebiet G,=G—G,, und erhalten analog zu Gl. (40) 


y/% 
(2.52) aytag f oF Gra lta; 0) O69; 7) oa = J galt xy) eo (A+ 9) dr 
G, 
mit 
(2.53) Gail”; 9) = ae [ (Go1 (3, £5 0) e, Xex) Vy |g — 0] 2K. 


ln—al=r 
SchlieBlich ist im Restgebiet G,: |3—y| 2y/2, d-h. 


Eee 
(2.54) | OF (3, 9; 27) e| =O(e : } 
daher . 
fe ae 
(2.55) amr | & Seles: 0) G3 (3,0; 47) e,dV, =O (e * ). 
G, 


Damit sind alle Ausdriicke abgeschatzt, die in der Bilinearentwicklung auftreten. 
Aus Gl. (50) bis (55) ergibt sich die Behauptung fiir Satz 2. 


§ 3. Anwendung eines Tauberschen Satzes 


Wir stellen zunachst einige bekannte Formeln iiber uneigentliche Integrale 
zusammen. Mit 


Y\x 


(3.1) Wt, %) =A = = Sinn ya — 2 f Ede 
gilt 0 
6-2) a ae dx=e"V({+y A) 
und 
(3.3) le(r, “| 2 (1+—) fiir alle 7, x=0. 
Ferner ist 
- 4 ia, 

ae ac tay ota wie V, 
und 

2 P sin s Va = 1 ite 
G5) n | ener i Yu ° rin tee O 


0 


Es sei r€G, ’,=a>0, und x eine nicht-negative reelle Zahl. Wir fiihren 
folgende Bezeichnungen ein: 


a/2 
6.) h(x,r) =JS g(r, x) u(r, x) dr, 
0 
und 
Gn (x)? x C P 
ae — —h(x,r 
G.7) A(*,t) =4 o<dns ye 0M me (x) —A(x,x) fir x>0 


0) fiir x=0. 
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Dabei sind C(r) und g(r, x) die bereits in Satz 1 auftretenden Funktionen. 


bilden wir 
4 f Ales) 
(3-8) B(A,x) af ee ax 
0 
und , 
C,, Z 
(3.9) Oty Se 
0<onSyx a 
Mit diesen Bezeichnungen gilt 
S1 IGP _ yf 2O(e2) 
een 1D) ot oh +A) Le: 
0 
bzw. nach partieller Integration 
< | En (x )P =e r 6, (%, x) 
(3.11) AO aoa i | aged 


0 
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Ferner 


Aus der Konvergenz der Stieltjes-Transformation Gl. (10) folgt naémlich bereits 


O, (x, x) =0(x), so daB die partielle Integration statthaft ist. 
Wir beweisen nun 


Lemma 3. Es 1st fiir x20 bzw. AZ0 und r,2a>0 


(3.12) Jha) =0(2), 

(3-13) | A(z, x)| =0 (x) +0 (—) 
und 

(3.14) |B(A,x)| =O(e 3"). 


Beweis. Die Abschatzung Gl. (12) folgt aus Satz1 und Gl. (3.3). 
ergibt sich nach Satz 1 und Gl. (3.11) 


(3.15) nee 


Da 0, (x, x) monoton wachst, ist 


(%, oH dx =0(\A) +0(-). 


lee) Cc 


(3.16) O,(2,2) =O, (22) 24 f 5 
A 


z+ A? = (v4)? 


Aus Gl. (7), (9) und (16) folgt Gl. (13). Nun ist nach Gl. (7) 


af i va) parma f 2 ee cs) mar 
yah i eleilcoeae 
Tia ofc a Grae” 


28* 


Ferner 


. 24 f n% 2) gy —O (A) +0(2). 
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Benutzen wir daher die Formeln (2), (4) und (6), so erhalten wir 


a/2 oo 


(3.18) ee Co) A fatead [eae axar 
aS 0 
= SR — Bow — felevene rr 
Nach Satz 1 ist mithin 
19 af Geeaenole 9), 
Andererseits ist ; : 
(3.20) A f ed dx =B(A,x). 


Aus Gl. (19) und (20) folgt Gl. (14). Damit sind sémtliche Behauptungen von 
Lemma 3 bewlesen. 

Weiterhin gilt 

Lemma 4. 


(3.21) ie *” A(v2, x) dv 


konvergiert fiir Re(z)>0 und ist dort holomorph. Die Funktion (z, x) ist holo- 
morph fortsetzbar in den Kreis |z|<a/3 und ist in |z|<a/4 beschrénkt. 


Beweis. Die erste Behauptung des Lemmas folgt unmittelbar aus Gl. (13). 


Ferner ist 
[ee 


\u 
(3.22) (et) = 5 _ " A(u,t) du. 


0 
Fiir reelle z>0 gilt daher unter Benutzung von Gl. (5) 


© 


cme f r Da esiieala” 
(3.23) pea=2f{ fator, Aaeue 4) ae. 
0 0 


Wir diirfen die Reihenfolge der Integrationen vertauschen und erhalten nach 
Gl. (20) 


2 Ree ) sing VA 
(3.24) p(t) =— f B(x) 22 aa. 
0 


Nun ist aber fiir komplexe z und reelle positive 4 
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d.h. 
|2| ya 7 a iD 
\ é a 1 —(Z—|z|) ya dA 
(3.26) je 2 [Bla dA o£ fe (%—lel) <). 


ya Va 


Die aoa rearsiftcts Gl. (21) und (24) stellen fiir positive reelle z dieselbe 
Funktion dar. Das Integral Gl. (24) konvergiert aber fiir alle |z|<«/3 und stellt 
dort folglich eine holomorphe Funktion (z, x) dar, die nach Gl. (26) in |z|<«/4 
beschrankt ist. Damit ist Lemma 4 bewiesen. 


0 


Wir benutzen nunmehr folgenden Tauberschen Satz von A. E. INGHAM}. 


Satz 3. Die Funktion A(x), A(0)=0, set reell fiir reelle x=>0 und von be- 
schrankter Variation in jedem endlichen Intervall. Ferner gebe es eine Funktion 
H(6)>0 so, daB fiir jedes 6>0 
(3.27) A(x'*) — A(x?) = — H(6) 


fiir xSx'Sx+6, x2 %(6)=0 ist. Das uneigentliche Integral 
(3.28) p*(2) = Je **d A(x?) 
0 


set konvergent fiir Re(z)=o>0 (z=o0+¢t), und *(z) erfiille die Bedingung: 
eo ish shetig-in o= 0,2 THIS eT. 
z 


Dann gilt 
(3.29) lim |A(x*) — 9*(0)| $M - A(T), 


B— iO 
wo M eine absolute Konstante ist. 
In unserem Fall ist 


A(x'2, x) — A(x?,r) = > [Gn (x)? ga {h(x'2, x) — h(x?, x)} 


oz is 


XSon< x’ 


(3.30) ue 
Te ae {u(r, %?) — u(r, x) } dr, 


IV 


und nach Gl. (3.1) 


vx 
sin € 
siny x’ — sinr x — {> zs as| 


Vx 


|m(r, 2) — u(r, 9) | = 


(3.34) 


so daB mit einer Konstanten 
a/2 


M,() = +4 [4-leealar 
pene aah de Se aponir Wece es, 3x0 
gilt. Ferner ist fiir Re(z) >0 
633) 9tle.t) = Jed Aa, 2) <2 fe AGA, 2) dv = 29,0). 


1 A. E. IncuHam [5], Theorem III dort. 
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Nach Lemma 4 ist also 


(3.34) aE SO 


stetig in co20, —a/4StS+a/4. Alle Voraussetzungen des Satzes von INGHAM 
sind mithin erfiillt und wir erhalten wegen y*(0, r)=0 als Ergebnis 


lim |A(x?,z)| SMM, ()-=. 


X—> CO (0 4 
Die Konstante M,(r) ist gleichmaBig beschrankt fiir alle r¢G. Es gilt daher 
Lemma 5. Es sev ¢€G, 7,20>0. Dann ist 


B (plz 2 
(3.35) O(ny= y MOP + sqny, 


(62) 
0<anSy% fe 


wobet S(x,x) der Abschitzung | S(x, x)| =O (1/u) geniigt?. 
Wir wollen aus Lemma 5 einige Folgerungen ziehen. Zundchst beweisen wir 
entsprechende Formeln fiir O;(x, x) (6==1). Es ist 


x 


> + |G, (x)? dO, (¢, x) 
3. 6 O x; r) = Ie. ee, 7 ? 
(3.36) olin Dh gs : ! 


und durch partielle Integration folgt 


x 


6.37) O (x2) = ED + (8-1) f C2 ae. 


Fihren wir fiir 0,(x, x) den in Lemma5 gefundenen Ausdruck ein, und setzen 
zur Abkiirzung noch 


S(t, x) Gert e—2) 
(6 yf ee (o>1, Es i) 
3.38) Wi) = = ; 
Dias dt 4 = (1 — log @,) (5=3) 
. 0 (<1) 


so erhalten wir folgendes Ergebnis: 


Satz 4. Es sei r€G, 7,24a>0. Dann ist fiir x—>0o 


pig i : : 
(3.39) O,(x, x) = saat el Olsen fiir b+ 5 
silogx +My) +0(12) fir 0=3., 


. — |G, (x)l2 3 . 5 
Insbesondere konvergiert Ss 1S5 (2) fiir 6 > = und divergiert fiir 6< ae 
et OND 2 2 


* Am Beispiel der Eigenschwingungen in einer Kugel la8t sich zeigen, daB dieses 
Ergebnis im allgemeinen nicht verbessert werden kann, — 
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Zum AbschluB dieses Paragraphen beweisen wir 


Satz 5. Die Anzahl der Eigenfrequenzen w,, im Intervall 0<w,<x geniigt 
fiir x> co der Abschatzung 


(3.40) | >» 1= —- + O(x? log x). 


0<anSx 


Huerbei ist V das Volumen des betrachteten Gebietes G. 


Beweis. Die Bilinearentwicklung Gl. (2.2) konvergiert fiir alle t=0, und 
zwar gleichmaBig im abgeschlossenen Bereich G. Setzen wir 


tim @2 (x, 2) =W,(x), 
so ist also W,(r) in G stetig. Andererseits ist nach Satz 4 


(3.41) 0, (x, 2) = Wa (e) — a + O(a). 


Oe A 


Wir integrieren diese Gleichung tiber G,, wobei wir zu G, alle die Punkte r€G 
zahlen, die von F einen Abstand 2@ haben: G,={r€G, r,20>0}. Dann ist 
zunachst 


(3.42) OS JO, (x7, 2) dV, S J Wa (et) dV, =O(0), 
G-Gy G-Go 
und es ergibt sich wegen der Stetigkeit von W,(r) in G 
, V 0 1 1 
“ae 0< ae [Mea Ax | Og) | o(%) | Ol; 108 >): 


Setzen wir jetzt 90=1/x?, so folgt fiir *— oo 


(3.44) Y BHM — a + 0 (log | 


n 
0<anS*x 


mit 


= [Wale ) aV,. 


Durch partielle Integration ergibt sich daraus 


(x? log x), 


_ V+ 
(3.45) l= tare 


0<anS% 


d.h. die Behauptung von Satz 5. 


§ 4. Weitere Ergebnisse 
Wir wollen in diesem Paragraphen zunichst beweisen, daB fiir y==y 


2 
On 


Ey (x) Ei (o) | =O(1) (x >) 


0<anS% 


ist. Dazu wenden wir den schon zitierten Tauberschen Satz von INGHAM auf 
die in Satz 2 gewonnene Abschatzung der Bilinearentwicklung an. Der Beweis- 
gang ist weitgehend analog zu den Ableitungen in §3. Allerdings werden wir 
an einer Stelle von dem Ergebnis des Satzes 4 Gebrauch machen. 
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Wir wollen annehmen, da die Punkte ry und 4 (rEG, HCG, r=y) fiir das 
Weitere fest gewahlt sind. AnschlieBend an Satz 2 fiihren wir folgende Bezeich- 
nungen ein (vgl. die analogen Bezeichnungen in §3): 


y/2 
(4.1) Ay 1 (%) = J built & b) u(r, x) ar, 
“hye OE 
0<onSy*% " 
(4.3) Axi (*) = O}%1(*) — Cyilt, 9) — Agi (), 
4 f Anl#) - 
(4.4) By, (A) =A tart | (k, 2 =4, 2, 3). 
0 
Dann gilt 
Lemma 6. Es ist fiir x0 bzw. AZO 
4.5 hy (x) 4), 
(4.6) Agi (x) = Of1(%) +0 (1) =O (x) +0 (1) 
(4.7) Byi(A) = 0 (62) 


Beweis. Die erste Abschatzung Gl. (5) folgt unmittelbar aus Satz 2 und 
Gl. (3.1). Gl. (6) folgt mittels der Schwarzschen Ungleichung aus Satz 4, denn 
es ist 


Jones y LAO. y Ix 


O<anSx On 0<anSy% on 

(4.8) < En (x) |? > |, (9) |? 
“aeaneye On 0<onSy% On 
= 0(x) 


Gl. (7) schlieBlich ergibt sich bei Verwendung von Satz 2 aus 


Bult) =A [ as dx 
; y/2 


ic) k 
a >) Za) Fal — Cy i(t, b) ae Eb) e- ed 6 +174) dr 


2 
. 
n=1 @ 


(4.9) 


=0 (e- 2 ie : 
Man vgl. hierzu die analogen Herleitungen in §3, Gl. (17) bis (20). 
Die Funktionen 


(4.10) Pri (2 afar *° A, (v?) dv (k,l =4, 2, 3) 


sind daher ebenfalls holomorph fiir Re (z)>0 und holomorph fortsetzbar in den 
Kreis |z|<-+/2 und beschrankt in |z|Sy/3. 


— 
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Wir wollen schlieBlich noch zeigen, daB8 die Funktionen A,,(x) auch die rest- 
lichen Voraussetzungen des Satzes von INGHAM (Satz 3) erfiillen und beweisen 
dazu 


Lemma 7. Es gibt positive Konstanten M,, und K, so daB fiir jedes d>0 


(4.11) Ay)(x'?) — Ay,(x*) 2 — (My, 0+ K) 


ast fiir alle x’ mtOSxSx’Sx+6. 
Beweis. Es ist 
, E& (y) El, . 
(4.12) Agila’) — Ay(xt) = SL Ee) Bal) _ fy (47) — ty, (29)}. 


, OF, 
X<WnS% 


Nach der Schwarzschen Ungleichung und nach Satz 4 bzw. Gl. (3.34) erhalten 
wir mit einer geeigneten positiven Konstanten K 


4<OnSx’ On at X< On Sx’ On Lay Se On 
(4.13) ; j 
(24% 
deh, a 
(4.14) hs En(t) En(o)_ aot ( mit K} 
wz = eee ‘ 


Ferner ist nach Gl. (3.31) und (4.1) 
y/2 
(4.15) hg s('®) — Ig (9) 2 — —— | - |Sar(es t,o) ar, 
0 
mithin insgesamt 


y/2 
Ay ;(x'*) — Ay, (%?) = —( 4. += [r|gurlrs x, y)| ir) aus 
(4.16) 0 


fiir OS xS<x’/<x+6 und alle 620. 
Die Voraussetzungen des Satzes von INGHAM sind daher erfillt und wir 
erhalten nach Satz 3 


(4.17) A, (x) =O) 
bzw. nach Gl. (4.3) und (4.6) 


0<anSx ie 


Wir wollen noch untersuchen, fiir welche Werte von 6 der Ausdruck 0%, (x) 


konvergiert fiir x—>oo. Es ist 
x 


d Oz (é 
(4.19) Oh) hy) = f SE, 


y 
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oder nach partieller Integration 


(4.20) O8,(x) — OF, (y) = 


Wegen O})(x) =O(1) ist 
: at = 1 1 
(4.21) |x, (x) — O::(y)| = o(—=) ts Oe , 
so daB die rechte Seite im Fall 6>1 beliebig klein gemacht werden kann fiir 
geniigend groBe x und y. @f,(x) konvergiert daher, falls 6>41 ist. Es gilt mithin 
Satz 6. Fir rEG, yEG, r=) und x00 1st 
Ent) Gi (y) =O(A1). 


oz 


0<a,S% 


Die Reihe >’ A honvergiert fiir 6>1. 
nai n 


Wir bemerken abschlieBend, daB wegen der Symmetrie der Greenschen 
Tensoren Gl. (1.12) aus der Bilinearrelation wie in Satz 4 und 6 auch fiir die 
magnetischen Eigenschwingungen gelten. Sind ¢ und w jedoch von 1 verschieden, 
so ergeben sich aus den Substitutionen Gl. (1.2) folgende Ergebnisse: 


Satz 7. Es seien &,, 9, die elektrischen bzw. magnetischen Eigenschwingungen 
zur Eigenfrequenz w, der Maxwellschen Gleichungen 


VxE—ZgouH=o0, Vxo+¢oeC=n 


mit der Randbedingung nx €=o auf F. Die Eigenschwingungen seien durch 
= [S, G,,aV =— [ Dn NS Onn 
G G 


normeert. 
Dann gilt fiir rCG, r,2a>0 und x00: 


> |16.@F= eit o(4 2) 


0<anSx% 
3 
2 
3 |, (x) = 228) x3 O( x). 
0<m,S%x % 
Die Reih IEn(2) 1m @) 
ie Reihen “ nd S  Ronvergieren fiir b> 3 - und divergieren 
wo; 
n=1 n=1 one 


fiir je 
2 
Ferner ist fiir y= und x->0o: 


Y 46.0 Gy) =011) 


ocores On 
ye Tats 
> 2 Dn (x) Dn ()) ara O (4) . 
eran 
BORE ~ 4 
Die Rethen >" “wns €,, (x) EF (y) und Dae 25 ,,(t) 52 (y) Ronvergieren fiir 5>1. 


n=1 
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Fiir die Anzahl der Eigenfrequenzen w, im Intervall 0<w,<x gilt die Ab- 
schatzung 


— 


i 1=(eu)3 — + O(x* log x). 


0<anSx 3 


Diese Arbeit entstand wahrend eines Aufenthaltes am Institute of Mathematical 
Sciences der New York University. 


[11] 


[12] 
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Electro-magneto-optical Effects 


R. A. TouPIN & R. S. RIVLIN 


1. Introduction 


Electro-magneto-optical effects in stationary materials are observed when 
the materials are placed in strong, static electric or magnetic fields and an electro- 
magnetic wave (light) traverses the medium. The Faraday effect is an example 
of the class of phenomena we have in mind. Clearly, the existence of electro- 
magneto-optical effects in material media is direct evidence that the equations 
for the electromagnetic field in a material medium, unlike the equations for it 
in vacuum, are non-linear since the sum of two solutions generally fails to be 
a solution. An analysis of electro-magneto-optical effects, though based on non- 
linear equations for which all but the simplest exact solutions are difficult to 
exhibit, may be simplified and made tractable by assuming that the dynamical 
part of the solution (the light wave) has such weak intensity that the vectors 
which describe it may be treated as infinitesimals. 

Accordingly, in this paper we assume that the electric displacement and 
magnetic intensity fields are functions of the static electric and magnetic induction 
fields and linear functionals of the small dynamical electric and magnetic induction 
fields. This starting point is derived, in the Appendix, from the assumption 
that for a more general class of fields, not necessarily representable as the super- 
position of small dynamic fields on large static fields, the electric displacement 
and magnetic intensity fields are non-linear functionals of the electric and 
magnetic induction fields. 

From our initial assumption regarding the variables which enter into the 
constitutive equations, we derive the limitations which are imposed on them 
by the assumption that the material is isotropic. The isotropic constitutive 
equations are then applied to the study of the manner in which a plane electro- 
magnetic wave is propagated in an isotropic material to which static electric and 
magnetic induction fields are simultaneously applied. 

From these results we consider in greater detail the cases when the applied 
electric field is zero, leading to the magneto-optical effects, and that when the 
applied magnetic induction field is zero, leading to the electro-optical effects. 

The case when both the static electric and magnetic induction fields are non- 
zero will be given more detailed consideration in a later paper. 


2. The constitutive equations for small varying fields 
superposed on large constant fields 
We assume that there exist in a body an electric field E* and a magnetic 
nduction field B*. Each of these is assumed to be the resultant of a time- 


0 a 
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independent field and a small field which may vary with time rv. Thus, 


E*(t) = + ¢ E(t) 
B*(t)=&+ ¢B(r), 


and (2.1) 


where & and & denote the time-independent fields and ¢ E(t) and é B(t) the 
small additional time-dependent fields. ¢ is a small constant. The electric dis- 
placement field D*(t) and the magnetic intensity field H*(rt) may also be ex- 
pressed in the forms 

D*(t) =D + e D(z) 


ee H*(t) =H 4 ¢H(c), en 


where D and # are the values of D*(r) and H*(t) corresponding to e=0. D 
and # are, of course, independent of time. For brevity we shall employ the 
notation 


D@)=D, H)=H, E()=E and Bi)=B, (2.3) 


where ¢ is a particular instant of time. 


In general, we may expect that for a given material the values of eD and 
eH will depend on @ and @& and on the values taken by eH(t) and ¢B(t) 
throughout the time interval —co<tSX?. We denote the components of D and 
H in a rectangular Cartesian coordinate system x by D,; and H; respectively, 
those of H(t) and B(t) by £;(t) and B,(t) respectively and so on. Also we use 
the notation 

i, (ie and. - B;(t) = .5;. (2.4) 
Then, each of the vector components ¢D; and eH; may be expressed as a func- 
tional of e£,(r) and eB, (t) over the range — co< 7S and an ordinary function 
of &, and B,, thus 


t 
sD, = Oleh, (2), eB, 6, Bp 
and rat (2.5) 
éH,— Vile E,(t), éB,(t); é,, By. 


— 


3. Restrictions on the constitutive equations 
for hemihedral isotropic materials 


If the material is hemihedral isotropic in the absence of any applied fields, 
then the constructive equations (2.5) must be form-invariant under the proper 
orthogonal group of transformations. Following a procedure similar to that 
adopted in previous papers [J, 2, 3], it can be shown that the restrictions imposed 
on ®, and ¥, by this form-invariance may be expressed in the forms 

t t 
O;[e Es (t), € By(t); &, By] = 4; Pile Ep (7), € By (7); $y, Hy] 
and cat ss (3.1) 
WieE,(t), e B,(t); &, B]=4,; Vile E(t), € By(t); &,,%], 
=—oo —oco 


t= Tt 
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where Eh a, Bvt) By ay Be) 
E> =4p,6, and By = Ay By, 
and a,, are real constants satisfying the equations 
Apo Ags= Asp sq =Opq, | Upq| =1- (3.3) 
We introduce the notation 
E, (t) =i B,(t)= Bye’, 
d* Ey (t) 


dt 


(3.4) 


—— fo) 


Pay eee ae (a= 1, 2,..., m4). 


In order to obtain the restrictions imposed on @®; and ¥ by the relations (3.1) 
in more explicit form, we assume that, in (2.5), ®; and Y% are continuous here- 
ditary functionals of «£,(r) and ¢B,(t) over the range —co<1Sz and poly- 
nomials in e£)” and eB) («=0, 1, 2,...,) andin & and Z,. We also assume 
that ¢ is small enough so that we may neglect terms in ®; and ¥Y% which are non- 
linear in e«. We use the notation 


B= (Dye) Meet aera) 
BO = (BM) for eres an eny 
Then, following a procedure similar to that adopted in previous papers and 
bearing in mind that ®@=W=0 when e=0, it can be shown that the vectors 


® and Y and hence the vectors ¢D and ¢B may be expressed as linear combi- 
nations of vectors of the forms 


and 


(3.5) 


t 
e EO, e f y(t — t) E(x) dx, 


eCXB®, c¢BxE, 


i (3.6) 
efx f y(t — 1) B(t) dt, 
t 
eBx f y3(¢ — +) E(t) dt; 
vectors of the forms oe 
t 
eB, és fy (¢— 7) Blt) dt, 
cOXE®, e&BxB, 
! (3.7) 
e€xX f ys (t — t) E(t) dt, 
t 
e Bx f x_(t — t) B(t) dt; 
é, »Ex® (3.8) 
and 
&, (3.9) 


where ce =A Ons Uy ere 7m 
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The coefficients of the vectors (3.6) and (3.7) are polynomials in 


f-E, B-B (3.10) 


and 


z- 8. (3.11) 


The coefficients of the vectors (3.8) and (3.9) are polynomials in (3.10) and 
(3.11) and are linear in the quantities 


Bo Bo 9 eB oBo, e[E™, €, B] (3.12) 
and 


ee: B®, ¢@-E@, -«[B,2,8], (3245) 
and in 2 obtained from Ue 12) and (3.13) by substituting vectors of 
the forms Sxt-2) E(t) dt and ja ¢—t) B(t) for the vectors E™ and B® 
Perec ae 

4. Restrictions on the constitutive equations 


for holohedral isotropic materials 


We note that the vectors (3.6) and (3.8) are changed in sign by the central 
inversion transformation while the vectors (3.7) and (3.9) are unchanged by this 
transformation. Similarly, the scalars (3.11) and (3.13) are changed in sign by 
the central inversion transformation, while the scalars (3.10) and (3.12) are 
unchanged. 

Since the constitutive equations for a holohedral isotropic material must be 
form-invariant under the central inversion transformation, in the constitutive 
equation for ¢D, each of the vectors (3.6) must have as its coefficient a polynomial 
in 

Gree Gerte and (€-76)?; (4.1) 

each of the vectors (3.7) must have as its coefficient an expression of the form 

at: d, (4.2) 

where « is a polynomial in the quantities (4.1); each of the vectors (3.8) must 
have as its coefficient an expression of the form 

eD{AME- EO + AMB- Be + AP EO, E, Bl} + 


+e€-BY{(CHL. BO + CH B- HO + CH (BO, &, B]}, (4-3) 


where the A’s and C’s are polynomials in &2-&, ®-& and (€-&)?, the vector 
(3.9) has as its coefficient an expression of the form 
BC Ha {FR &-BO 4 FOB. Bo + FOOTE, , Bi ze 


sexs Sy {GW é. Bo Ap Ge) bB. BO + Ge (Bo, & Bj, (4.4) 


where the F’s and G’s are polynomials in &-&, ®- & and (€- &)?. 


Again, in the constitutive equation for eH, each of the vectors (3.6) must 
have as its coefficient an expression of the form (4.2); each of the vectors (3.7) 
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must have as its coefficient a polynomial in the quantities (4.1); each of the 
vectors (3.8) must have as its coefficient an expression of the form (4.4) and the 
vector (3.9) must have as its coefficient an expression of the form (4.3). 


5. Constitutive equations for sinusoidal time-dependent fields 
We now consider that the fields ¢H(r) and eB(r) vary sinusoidally with 
time with angular frequency w, so that we may write 


E(t) =2E(x)e'®* and B(t) =Z B(x) ec", (5.1) 


where a is the position vector of a generic point of space. We shall, for the 
remainder of this paper, follow the common convention by omitting the symbol &. 
Then, introducing (5.1) into (2.1), we obtain 


E*(1) =2 + e E(x) e'°* 


B*(t) =B+ ¢ B(x) e~'*. ee 


and 


Introducing (5.1) into the constitutive equations for a holohedral isotropic material 
discussed in §4, we obtain 


D(t)=D(x)e°* and H(t) =H (ax) e~!**, (5.3) 

where D and H are given by expressions of the form 

D=H4,E+a8XE+%€xXB+ 
+€-B (a,B+a,€XE+ 0,8 xB) + 
+ {a,E-E+a,8-B+ mE, , BE + 
+ {ag E+ op BB+ cry) [H, €, BY ExB + 
+2€-B lay E+ Bt oy. B-E+ a1 [B, &, BY} E+ 
+E-B{o,€- B+ 4B-E+ tt. [B,&, B}ExB+ 
b {y5E+ B+ a6 B- E+ 0% (B, f, Bl} B+ 
+E Blo + E+ a4, 8- B+ ay, [H, €, BB 
H=6,B+p,€xE+p,8xB+ 
+&-B(6,E+6;8xE+B,exB) + 
+{6,€- B+ py B-E+ Byo[B, €, BE + 
+ {65 E- B+ By B-E + Boo [B, 2&2, BExB + 
++ B {By E-E+ BBB+ By [E, ef, BY E+ 
+E BiB é- E+ By B- B+ By (B,E, BYE x B+ 
+ {Bis€- E+ Big ®- B+ Bos (E, é, BB+ 
+é- Bib, €- B+ Bis 8-E+ B4(B, 2, BB, 


Where the «’s and f’s are functions of ww and polynomials in &- f, @8-& and 


(C- B)?. 


| 


and 
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The following identities hold: 
[E,e,B\E=—(€-BExE+ (€-OBxE+ (C-BE)Ex&, 
[B, ¢, B)€=—(€. 8) ExB—(E-2)BxB+ (€-B)ex&, 


[E, 0,8) 8 = (£-8) BxE—(&-BExE— (B-BE)Ex&, Oe) 
[B,€,B)B8= (6-8) BxB+(&’-8)ExB-(B-B)ExB. 
From the third of these we obtain 
[E,¢,B)ExB = (€-&B)Ex(BxE) — (8-B)Ex(ExE) — 
—(B-B)Ex(Ex&B). oe 
Using identities of the type 
&x(B XE) = (€- E) B—(é-&)E, (5.8) 


we obtain from (5.7), 


[E,2,8)ExB = (é-&) ((€-E) B— (€-&) E] — 


— (8-&)[(- £)& — (€-&) E] — (5.9) 
= (&- Ey((é-@)E—(€2) &]. 
Similarly, we have 
[B, 2, B)ExB = (~- 8B) ((€- B) &—(€- &) B) — 
= eb (Cab) C— (cc) Bl — (5.10) 


=(EoBy iC 6) C— (Co) 6). 


Using the relations (5.6), (5.9) and (5.10), it is clear that we may, without 
HissuOle Cenerauty, take” 64, Gg —"" =o, — 0) aNG® Py, —=Pa9 = = Pag Om 
equations (5.4) and (5.5). 


6. Constitutive equations for a plane wave 


We now consider that the time-dependent fields are those for an infinite 
plane wave, with angular frequency w and wave number k, in the direction of 
the unit vector m, and are given by 


Bie en em?) 
B=belltna-on), 
D=dein-a—o8), 
H=hetr ead, 


where the complex vectors e, b, d and h are independent of position and time. 
a denotes the position vector of a generic point of space. 


In this case, we may replace E, B, D and H by e, b, d and h in the con- 


(6.1) 


stitutive equations (6.1), in which we take 049 =o) =-+:=%4=0 and Bg =Bo9 = 
---==B5,=0. These equations may then be re-written as 
d=@.e+W.-b, 
x (6.2) 
h=Q-e+A-b, 
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re 


where ®, W, Q and A are the 3 x3 matrices defined by ®=|@;,;|, P=|¥%; 


and 
D, , = 0 0; 3 — Xp Erg By — Ms Em Bor, Fi je Ep + Oy G; G + 
+ XK Eimn &m Bn Oj + M2 Em Bm €;, B+ 


+ 034 En Bm Fin &p B, B + Me BB+ 
+ Oy En, Bn B; &; 


Yj = — Og €;54 SG + ha En By 0:5 — 6 Em Bry Exp By 
+ ag &; B+ O10 Eimn Sm By B + 
+ O41 En Br ©; & + O43 En Bn Eipg Sp B Fj + 
+ 045 B; & + 053 6, Bn BF; (63) 


i; = — Bo bizn E+ Ba En Bn 9:5 — Bs Om Bn Fig Be + 
+ Bs €; F + Bro Eimn Em B, Bj + Bu Ey By 6; ; + 
+ By; Em Bn Eipg é, B, €; + Bis BE, + 
+ Bis $n, By, B; B, 


A;; = py 6;; — Bs Eijk B, — Be Em Bn Eijk é;, + B, é; 6; + 
+ By Eimn Em B, €; + Bis En By E, B+ 
+ By En B,, Epa E, 2, B+ Big BB + Bir En B,, B, E,. 


6;; denotes the Kronecker delta and ¢;;, the alternating symbol, and the a’s 
and f’s are functions of 1 and polynomials in &-&, & -® and (€- &)?. 


7. The propagation of a plane electromagnetic wave 


The electromagnetic field equations in the absence of free currents and 
charges are 
0.B* : 
37h Corl BY 0} dive = 05 


(7.0) 
curl H* — FX div D* =o. 


We now assume that E* and B* are given by (2.1) and that the static portions 
of these fields, & and @&, are independent of position. We then see that FE, B, D 
and H must satisfy the equations 


St cuwlE=0, divB=0, 


(7.2) 
curl H — £2 =o, div D=0. 


For the infinite plane wave (6.1), we obtain 


knxe—wb=0 and knxh+od=0. (7.3) 
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Introducing (6.2) into the second of equations (7.3), we obtain, with the 
notation * 7 =k/w, 
nnx(Q-e+A-b)4+@-e+W¥-b=0. (7.4) 
Using the first of equations (7.3) to substitute for b in (7.4), we obtain 
yn x[Q-e+nA.- (nxe)]+nW¥-(nxe)+@-e=0. (75) 


(7.5) is a system of three homogeneous linear equations in the components of e. 
We may therefore write it as 


177 =9; (7.6) 
where 
RPS @;; Sas) (Ej x1 Yn — E;p1 2p; m) +? EjkI Eipg Lg Mp N- (7.7) 
The secular equation for 7 is 
[xij] =0. (7.8) 


If the wave is propagated in the direction of the x-axis of a rectangular 
Cartesian coordinate system x, we may take 1;=63; in (7.7). We then obtain 
from (7.8), 


G49 (Ae— 221) — 7? Age, Dy o— 7 (Yt Qo) +7? Ap1, D,3— 7? Qo 
Doi +7 (YetQ211) +7? Aye, Doo+7 (Yi t21.) — 7? An, Dy3+ Qi3 =0. (7.9) 
©, +7 ¥%., Dz. — V1, Ps 


This is the secular equation for an infinite plane electromagnetic wave pro- 
pagated in a holohedral isotropic material to which static electric and magnetic 
induction fields are applied. We note that it is an equation of fourth degree 
in 7. In general the solutions of this equation are complex. 


Denoting 4 by Ha ae. 7.10) 
where 7* and 7 are real, we see that a solution with 7*>0 corresponds to a 
wave travelling in the positive direction of the unit vector nm, while a solution 
with 7*<0 corresponds to a wave travelling in the negative direction of n. 
The amplitude of the wave will increase in the direction of propagation of the 
wave for the first of these waves if 7-<0O and for the second of the waves if 
n->0. Consequently, if either of these situations arises for any value of a, 
the system will be unstable. We shall in our subsequent discussion assume 
that such situations do not arise. 

Thus, the four solutions of (7.9) correspond, in general, to four elliptically- 
polarized, skew E-waves travelling with different complex velocities. The manner 
in which the plane of polarization and ellipticity can be calculated is shown in 
the next section. 


8. Analysis of a skew, elliptically-polarized plane wave 
We consider a plane wave given by 
E=ee'?, (8.1) 


where g=kn-x—ot. (8.2) 


* 1/n is then the complex velocity. 
29* 
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a 


We have* E*=e* cosy —e' sing. (8.3) 


The vector Et therefore executes an ellipse. Let I be a unit vector which 
is at all times perpendicular to the plane of the ellipse. J is taken in the same 
sense in relation to the plane of polarization as n, 7.e.1-n20. Then, l- E*=0 
From (8.3), we then have 


Lietca lt Qrand, (8.4) 
whence 
ike — * 
ip CK e | — exe (8.5) 
" (et x expaler xen j}i4 - [(exe*)- (ex e*)] 


the positive or negative signs being chosen to make l-n20. From (8.3), we have 
E*. E* = e*- e+ cos? y — 2e*-e-cosysing +e -e sin’. (8.6) 


We note that E*- E* is the square of the radius vector to a point on the ellipse. 
The vector E* is in the direction of an axis of the ellipse when d(E*- E*)/dg=0O, 


1.e. when 


2e*-e- t(e- e—e* - e* 3 
tan 29 = —; =o — an : (8.7) 
e*-e*—e -e- e-e+e*-e 


From (8.6) and (8.7), we see that the semi-axes L, and L, of the ellipse are given by 


TA, Test (Oe op eee) Se IG 3G a eines ae 
1 5 
2 


e-e* + 1f(e-e) (e*- e*) |). 


l and e may satisfy certain special relations which correspond to the KE*- 
wave having a particular character. These relations together with the nature 
of the wave to which they correspond are listed below: 


[e*,e,U]>0 right-handed wave, (8.9) 
[e*,e,I]<0 _ left-handed wave, (8.10) 
e-e=0 _ circularly-polarized wave, (8.11) 
e*xe=0 _ linearly-polarized wave, (8.12) 
t=n, andhence e-n=0O, transverse wave, (8.13) 
exn=0 _ longitudinal wave. (8.14) 


9. Magneto-optical effects in holohedral isotropic materials 
In this section, we shall discuss the interaction between a plane electro- 
magnetic wave in a holohedral isotropic material and a static magnetic field 
applied to the material. We accordingly take €=0. From equations (6.3), we 
then obtain 


D;, = % 0;; — Hy & 3, B + %5 FB Z,, 
A; = By 0;; — Bs &:j% B + Bre BF; 


* Here and subsequently, we shall denote the real and imaginary parts of a 
complex quantity by using the superscripts + and — 


(9.1) 
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and Yj, =2;,=0, 


v 


the «’s and f’s being functions of uw» and polynomials in &, Z,,. 
From (7.7) and (9.1), we obtain 


Mig = Diy EN? er ei pp Age Mp (9.2) 
a) Longitudinal magnetic field 
If the magnetic induction field @ is parallel to the direction of propagation 


of the wave, then we may without loss of generality take both of these directions 
as the 3-direction of a rectangular Cartesian system x. We then have 


With (9.3), the expressions (9.1) for ®;; and A;; become 
D;, = C4 07; —% B Exj3 + 6 B? 0s ; 03; 
A;; = B, 0; —p,B E573 + Big F? 0g; Os js 
where the «’s and #’s are polynomials in 4. From (9.2) and (9.4), we obtain 
Nig — %, Oi; — %B E4j3 + Ge B? dz; 03; Spe jolla on (0; — 03; 03 ;) + p3B E73] . (9.5) 
Whence, 


and (9.4) 


oo PH, <= (62 — P37), 0 
xii =||(% — B37?) Z, a, — B77, 0 . (9.6) 
| 0, 0, oy + Hy B 


The secular equation (7.8) then yields 
+ 04 B= 
Ee OT 46 (9.7) 
(o%, — Bi”)? + (o2 — B37)? B? = 0. 


Discounting the first alternative, we obtain 


_ atta, B 1 —1 by B (9.8) 


2 
1 BAR "Bip a 
From (9.6) and (7.6), it is easily seen that es=0. Consequently, since 1;= 0s ;, 
the relation (8.13) is satisfied. Also, unless (% ,—f3n?)#=0 or «,—f,n?=0, it 
is seen that the relation (8.11) is satisfied. Thus each of the four solutions for , 
given by (9.8), corresponds to a transverse, circularly-polarized plane wave. Two 
of these travel with speeds, which are in general different, in the positive direction 
of the x-axis and the remaining two travel with these speeds in the negative 
direction. If #>0, the wave in the positive direction of the x-axis given by the 
first of the solutions (9.8) is a left-handed wave and that given by the second 
of the solutions (9.8) is a right-handed wave; on the other hand the waves in 
the negative direction of the x,-axis given by the first and second of the solutions 
(9.8) are right-handed and left-handed respectively. 
It is thus seen that the right and left circularly polarized E* waves, with a 
given direction of propagation parallel to the magnetic induction field &, travel 
with different speeds. This is equivalent to showing that the plane of polarization 
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of a linearly polarized wave rotates at a rate proportional to the distance travelled. 
This is the well-known Faraday effect. In (9.8), «, B,, %: and fs, are, in general, 
polynomials in B?. 
b) Transverse magnetic field 

We now consider that the magnetic induction field @ is perpendicular to 
the direction of propagation n of the plane electromagnetic wave. We choose 
a rectangular Cartesian reference system with its 1-axis parallel to @& and its 
3-axis parallel to m. Then, in (9.1) and (9.2), we take 


B=Bs,, and n,;=0;;. (9.9) 
We obtain 
Dig = 03. Oj —= Git Spek ip reap: (9.10) 
A;; = B1 5;; — Bs F &i;1 + Big F? 61; O1;, 
and 
ij = %1 05; — hp B &;51 + yg B? 0, ; 04; + (0.411) 
+? [— B,(0;; — 63; 03;) — Big B? &:13 113] - 
Whence, 
Oy + O16 B? — By 7?, 0, 6) 
lal= 0, a — 4°(B, + Bis B?), —o%Bl. (9.12) 
0, ty B, Oy 
The secular equation (7.8) then yields 
2 % TK. B A 
7 ” 3. (9.13) 
an Lie 
EL ea 14 
4 4 (B, + By 6 B?) 0 ) 


If we take the solution (9.13) for 7?, we see from (9.11) and (7.6) that, in 
general, ¢,=e,;=0. Thus the two waves corresponding to this solution satisfy 
the conditions (8.12) and (8.13) and are therefore transverse waves, linearly- 
polarized in the x,-direction, travelling with equal speeds in the positive and 
negative directions of the x-axis. 

The solution (9.14) for 7? yields, with (9.11) and (7.6), e,=0. It therefore 
corresponds, in general, to two elliptically-polarized skew waves, the planes of 


polarization of which are normal to @&, travelling with equal and opposite 
velocities. For these waves, 


€3[€, = — & Blau, (9.45) 
so that e has the form 
e=£(0, 1,¢); (9.16) 
where 
C= — Oy Blory (9.17) 


and € is an arbitrary complex number. The Et-waves whose velocities are 
determined by (9.14) are linearly-polarized if and only if e* xe=0 (cf. equation 
(8.12)). This condition is satisfied if and only if € is real. The waves are trans- 
verse if and only if /;=0 (cf. equation (8.13)), i.e. if and only if €=0. They are 
circularly-polarized if and only if e.e=0 (cf. equation (8.11)), de. if and only 
if €=-tu. For all other values of ¢ corresponding to the solution (9.14) for ne 
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the E*-wave is an elliptically-polarized skew wave. The locus of E* at a fixed 
point in the material is an ellipse in the plane normal to &. From (8.8) the 
semi-major and semi-minor axes L, and L, of this ellipse are given by 


Li, L3=4E {1+ O0* + [(1 +02) (1 +0*2) #3, (9.18) 


the positive sign applying to the semi-major axis L, and the negative sign to 
the semi-minor axis L,. 


Introducing (9.16) into equation (8.7), we obtain 
nce dares eas taliao tT): 
ed TE od eh sete a ide) 
VoictT [5] explained the double refraction of light propagating normal to a 
magnetic field by the difference between the speeds as determined by (9.13) 
and (9.14) and observed the effect in vapours. Corron and Mouron (see, for 
example [6]) observed a very strong effect of this nature in certain liquids. 
The skewness and ellipticity of the E*-wave propagated with speed given by 
(9.14) seems not to have been noticed. 


tan29 = (9.19) 


10. Electro-optical effects in holohedral isotropic materials 


We now discuss the interaction between a plane electro-magnetic wave and 
a static electric field in the material through which the wave is propagated. 
We accordingly take ®=0. From (6.3), we then obtain 


D;; = 0 6,5 + oy 6; 6; 


B= — OX EjiR Ex» 


(10.1) 
Q,; =— B, EGR oe 
A;; =B,0;; +B, & &, 
the «’s and f’s being functions of ww and polynomials in &, &,. 
From (7.7) and (10.1), we obtain 
Lig = (01 O53 + Hy &; &) — n[(Bo + %) Ey, Mm Oi7 — (Hg 6; + By &;7;)] + (10.2) 


+ 7? [Bi (1; OL 0:4) + By EGR Ei pg 6, Ey Mp nm]. 


a) Longitudinal electric field 
If the electric field & is parallel to the direction of propagation of the wave, 
then we take this common direction as the %,-axis of the coordinate system %. 


Introducing 
n,;=6;, and 6,=€ 063; (10.3) 
into (10.2), we obtain 


Lig = (Ota 6;; + Oty &? 65; 63;) — €(Be + %) (5;; — 931 53) — 9° Br (0:3 — 43; 53 ;)- (10.4) 
Whence, 
oy — (By + %) &n — Bi 7?, 0, 0 


las, = 0, Ora (Bo + Gg) En a Bin? 0 : (10.5) 
0, 0, &, + a, & 
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Discounting the possibility that «,-+-«,6=0, the secular equation (7.8) then 


ield 
yields — &(B,+%) +[ (By +o)? 62+4 m% pt (10.6) 


2B; 

Unless B,-+-4,;=0, these two solutions for 7 correspond to waves travelling with 
different complex velocities parallel to the xs-axis. We may expect that one 
of these waves will travel in the positive direction of the x-axis and the other 
in the negative direction. 

Introducing (10.5) into (7.6), we obtain e,=0. Consequently, for each of 
the waves the relation (8.13) is satisfied and the wave is accordingly transverse. 
Also, it is easily seen that the solutions (10.6) for 7 give the inverse complex 
velocities for waves linearly-polarized in an igvaenin direction normal to the 
direction of propagation. 

We note, as was pointed out in §7, that if 7*>0 and 7 <0, or if 7*<0 
and 7~>0, the system is unstable. Consequently, %,, £,, B, and «, must be 
such that this situation does not arise for either wave. A number of conclusions 
regarding the values of «,, 8,, Bg and «,; might be drawn from this fact by purely 
algebraic reasoning. For example, the situation cannot exist in which [ (B+ 3)? 
+ 4a, 8, ]3/B, is real and &(B,+«,)/8, is purely imaginary. 


Uae 


b) Transverse electric field 


We may, without loss of generality, choose the reference system with the 
%,-axis parallel to the direction of wave propagation and the x-axis parallel 
to & We then write 

n;,= 03; “and” €,== ¢ 6,;- (10.7) 


Introducing these relations into (10.2), we obtain 


Nig = (Oy 05; + oy CP by; Oo;) +77 & (%3 0g ; Oo; + Be Oo; 63;) + 


(10.8 
+ 7? [B; (03; 03; — 0;;) + By € (do; 0»; + 63; 0s; a 6;;)] . ) 
Thus, 
— 4? (By + By 6), 0, 0 
xl = 0, O% +t 6° — Bin”, BE (10.9) 
0, Xs én, Oy 
The solutions for 7? of the secular equation (7.8) are therefore 
2— 
BB; &’ eg 
or 
ese Oy (% + Oy &*) (10.11) 


ty B, + a5 Bo & ; 


The solutions (10.10) for 7 correspond to waves travelling with equal speeds 
in the positive and negative directions of the x,-axis. Unless 


Oy By 6? = x3 By — By aq — a4 By (10.42) 


€2—=€g3—=0 for each of these waves and consequently they are transverse waves, 
linearly polarized in the x,-direction. 
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The solutions (10.11) for 7 also correspond to waves travelling with equal 
speeds in the positive and negative directions of the x3-axis. For these waves, 
unless the relation (10.12) is satisfied, it is seen from (7.6) and (10.9) that 


0 
and : 


(10.13) 


Cg[€y = — Os E n/t. 


The waves are consequently skew and elliptically polarized in the x x5-plane. 
From (10.13), we see that e has the form 


e=£(0,4¢), (10.14) 
C=— % 6 nloy (10.15) 
and € is an arbitrary complex number. Following the discussion in §9, we again 
see that the semi-axes L, and L, of the ellipse executed by the vector E* are 
given by (9.18) and the values of the parameter w determining their direction 
are given by (9.19), where € is now given by (10.15). 


where 


11. Appendix 


We suppose that the magnetic intensity and electric displacement fields, 
H* and D*, at time ¢ are functionals of the electric and magnetic induction 
fields, E*(rt) and B*(r), over all times preceding the instant ¢. Thus, 

t 
D* = ®[E*(r), B*(r)] 
and ae (11.1) 
He PE (a), Ba) |. 


We assume that the functionals ® and Y satisfy the relations 


@[E*(t — i), B*(t —h)|=® [E*(z), B¥(1)] 
eee Seer (11.2) 
WlE* (7 — h), BY —h)| EBT) 


The material considered then possesses invariable heredity. 


and 


We also assume that each of the functionals @ and Y may be particularly 
dependent on the values of the arguments and of their time derivatives of various 
orders at time t=? and is continuous for —co<t<t?. Using the notation 

E*(t) = E*¥ = E*), Et — [d" B*(z)/a tv") 3, 
Bi(i)\= Bt — B* 0) Bt (4) — [a B*(t)/dt"),-1, (11.3) 


(ud 2, es) 


we suppose that the functionals ® and Y may be expressed in the forms 
ty t 
®@[E* (1), B*(t)] = 2 f - FLE*(1), B*(x)], 
fe Ce (11.4) 
Y[E*(r), B*(t)] = 2g G[E*(1), B*(z)], 


T= — CO T=—C 


where f and g are differentiable functions of E*™ and B*) (w=0, 1, 2,...) and 
F and G are continuous functionals of their arguments. 
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We now assume that the fields E*(r) and B*(r) are the resultants of static 
fields & and & and small varying fields eH (t) and eB(r), where ¢ is a small 
constant, as given by equations (2.1). We have, neglecting terms of higher 
degree than the first in e, 

i . é af : 
{(E*™, B*)—f(E*, Be )|--ote ys { om aaeor i BM | 
“u=0 | 


and va (41.5) 


* 


E+ 
0 


e= 


t t nt 
FE*(z), B*(t)] =F(@, 8) +e f p¢—1) E(t) dt +ef p¢—1) Blt) dt, 
where @ and w are functions of & and &. Here, we use the notation E=—E, 
BO —B, E™ =[d“E(1)/dt"),-;, BY =[d“ B(t)/dv"),.,. We also use the notation 


t 
F(¢, 8) = FE*(2), B*(7)]|.n0, (11.6) 
observing that &(=E*(t)|,.9) and &(=B*(r)|,-») are time-independent fields. 
The second of the results (11.5) follows directly from the integral representation 
of a continuous functional. We note that df/0E*™|,.) and df/éB*™|,_» are 
functions of & and &. 
It follows from (11.4) and (11.5) that, neglecting terms of higher degree than 
the first in e, the functional ®[H*(rt), B*(t)| may be expressed in the form 
t “ t 

@[E*(r), B*(t)] = O(E, B) +e) (a, E™ + 8B, B”) + f p(¢—1) E(t) dt + 
mas ay oie (11.7) 
+ fp ¢é—1) Blt) dt, 


where the a’s, B’s, p and w are functions of & and & which are, in general, 
different from those in (11.5). 
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